AD- A 198  776 


■  -  _  ■  <  -  F  - 

REPORT  DOCUMENTATION  PAGE 

- ~=r  ,  ..  ,4. 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

I.  REPORT  NUMBER 

AFIT/CI/NR  88-  1^1 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  TITLE  fandSubllllcJ  CH  A  ft  £.tL  1  A  Ol^  <3clt\S  j 

)4tT£ft£>37(l.uc70fLtS  dy  vA(l/A8L6  AuGif 
5p  ecTftO  SCO  piC.  ELL  1  p5o  M^Tfty 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

MS  THESIS 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AuTHORfJi; 

K^^0^3£7H  G.  MtftKCL  IL 

B.  CONTRACT  OR  GRANT  NUMBERf*; 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

AFIT  STUDENT  AT:  0 Ul  M Hft 5 ; Ty  OA  lOieilPiSK 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREm  a  work  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  report  DATE 

1988 

13.  NUMBER  OF  PAGES 

133 

4.  MONITORING  AGENCY  NAME  4  AOORESSf/l  d///eren(  from  Controlllnt  Oltlce) 

AFIT/NR 

Wright-Patterson  AFB  OH  45433-6583 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

I5«.  DECL  ASSIFICATION/DOWNGRADING 
SCHEDULE  ^1^ 

_ _ k _ 

6.  distribution  STATEMENT  ("o/ R«PO*‘0  <■.  ^ 

DISTRIBUTED  UNLIMITED:  APPROVED  FOR  PUBLIC  RELEASE 

w 


17.  distribution  statement  (of  thv  abattaet  entered  In  Block  20,  It  dlllereni  from  Report) 

SAME  AS  REPORT 


18.  SUPPLEMENTARY  NOTES 


Approved  for  Publ  ij^Releasei  lAW  APR  190-1 
LYNN  E.  WOLAVER 

Dean  for  ResearCT0md  Professional  Development  ' 

Air  Force  Institute  of  Technology 

_  Wriqht-Patterson  APB  OH  45433-6SB3 _ 

[l9.  KEY  WORDS  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 


I  20.  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

ATTACHED 


88  9  6  20  e 


DD  ,^°r73  1  473  EDITION  OF  I  NOV  65  IS  OBSOLETE 


L'rJCLA_SSOT 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  flWien  Date  Entered) 


The  Characterization  of  GaAs/AlGaAs  Heterostructures 
by  Variable  Angle  spectroscopic  Ellipsometry 

Kenneth  G.  Merkel  II,  M.S. 

University  of  Nebraska  -  Lincoln,  1988 

Advisor:  John  A.  Woollam 
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thickness,  alloy  composition,  and  growth  quality  of 
GaAs/AlGaAs  samples.  The  following  separate  studies  were 
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roughness  measurement.  All  samples  were  grown  by  molecular 
b^am  epitaxy  (MBE)  at  other  research  facilities. 

^  In  the  first  study  the  layer  thicknesses  and /'alloy 
compositions  of  a  sample  were  measured  on  a  multilayered 
sample  at  three  positions  extending  radially  from  the  cen-  ^ 
ter  of  the  sample.  Thus  a  three-dimensional  analysis  of  the 
sample  was  performed  over  a  one  inch  area.  The  layer  thick¬ 
nesses  of  the  AlGaAs,  GaAs,  and  native  oxide  were  found  to 
vary  by  less  than  one  percent.  No  evidence  of  an  interfa¬ 
cial  roughness  layer  was  detected  on  top  of  the  AlGaAs. 
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I.  Introduction 


Rapid  progress  has  occurred  in  the  development  of 
GaAs/AlGaAs  modulation-doped  and  multiple  quantum  well 
heterostructure  devices.  These  types  of  structures  have 
found  applications  in  opto-electronics,  digital  electron¬ 
ics  and  microwave  analog  communications  (1). 

These  advances  have  created  a  requirement  for 
monitoring  the  accuracy  and  quality  of  III-V  semiconductor 
growth  processes.  Variable  angle  of  incidence  spectro¬ 
scopic  ellipsometry  (VASE)  provides  an  accurate  and 
nondestructive  method  of  characterizing  both  simple  and 
complex  multilayered  semiconductor  structures  (2) .  Par¬ 
ticularly,  VASE  is  applicable  as  a  diagnostic  tool  for 
studying;  the  layer  thicknesses  and  alloy  compositions 
(3),  wafer  homogeneity  (4),  surface  oxide  growth  (5),  in¬ 
terfacial  roughness  (6),  and  electric  fields  (7)  in 
heterostructures  and  superlattices. 

Chapter  1  contains  an  overview  of  the  VASE  process 
including  the  theory,  measuring  system,  and  data  analysis. 
Chapter  2  provides  a  background  discussion  of  the  physics 
and  technology  of  GaAs/AlGaAs  heterostructures  in  areas 
pertinent  to  this  research  effort.  Chapter  3  encompasses 
the  results  of  VASE  studies  performed  on  GaAs/AlGaAs 
samples . 

A  note  to  the  readers:  each  section  of  this  thesis  is 
a  self  contained  unit.  References  are  provided  at  the  end 
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of  each  of  these  sections.  Figures,  tables,  equations  and 
their  variables,  unless  explicitly  mentioned,  pertain  only 
to  the  section  in  which  they  are  mentioned. 
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Chapter  1. 

variable  Angle  Spectroscopic  Ellipsometry 


A.  The  Theory  of  Ellipsometric  Measurement 


The  theoretical  analysis  of  an  ellipsometer  system 
containing  a  polarizer,  sample,  and  analyzer  (PSA  system) 
is  fairly  simple  in  concept.  The  VASE  apparatus  can  be 
represented  in  a  manner  analogous  to  any  physical  system. 
A  "transfer  function"  for  the  VASE  is  derivable  in  terms 


of  the  effect  of  each  optical  component  on  the  input  elec¬ 
tric  field.  An  elegant  and  utilitarian  approach  to  any 
ellipsometric  configuration  is  available  through  the  Jones 
vector  and  Jones  matrix  formalism  (1) 

A  representation  of  the  amplitude  and  orientation  of 
an  electric  field  is  provided  by  the  Jones  vector.  The 
initial  step  in  the  Jones  approach  is  to  place  the  scalar 
components  of  the  electric  field  in  a  vector.  These  sca¬ 
lar  components  can  be  expressed  in  terms  cartesian  coordi¬ 
nates  by  p  —  “n  P" 

♦  r>))  (1) 


This  is  an  electric  field  propagating  in  the  +z  direction 

with  component  amplitudes  E  and  E  ,  angular  frequency  w, 

X  y 

wavevector  k,  and  phases  6  and  6  .  The  temporal  aspects 

X  y 

of  the  wave  can  be  suppressed  since  the  field  components 
for  a  monochromatic  wave  are  known  to  oscillate  in  time  at 


the  same  frequency. 
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By  considering  the  field  with  respect  to  a  single  plane 
perpendicular  to  the  direction  of  travel  (i.e.  a  plane 
situated  at  z  =  0) ,  the  spatial  dependence  may  also  be 
dropped  -  Te.  “I 


r" 

I 

i 


e  ' 

£u 


This  is  the  Jones  vector  representation  of  the  electric 
field. 

The  Jones  matrix  describes  the  overall  effect  of  an 
optical  system  (or  the  system's  components)  on  an  incident 
wave.  A  schematic  diagram  of  an  optical  system  with  a 
plane  wave  input  for  cartesian  coordinates  is  shown  in 
Fig.  1.  The  input  wave  has  the  right  handed  coordinates 
x,y  and  z.  The  output  wave  has  coordinates  x',  y'  and 
z'.  The  directions  of  travel  are  z  and  z'  for  the  incoming 
and  outgoing  wave  respectively.  S  signifies  the  optical 
system  (i.e.  the  sample  in  the  PSA  system).  In  the  absence 
of  nonlinearities  and  frequency  altering  effects,  the  out¬ 
put  wave  components  in  the  (x',  y',  z')  reference  frame 

can  be  designated  in  terms  of  the  input  wave  components  in 
the  (x,y,z)  frame  in  the  following  manner: 

-  T\i  *♦*  Til 

-  Tj,  £([;(  122  ECjji 

These  two  equations  are  easily  put  into  matrix  form 


rGo*'!. fT,, 

u"'«-  -1 


Note  that  equation  (5)  gives  the  Jones  vector  for  the  out- 
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put  wave  with  respect  to  the  Jones  vector  of  the  input 
wave.  The  transition  matrix  between  the  input  and  output 
Jones  vectors  is  the  Jones  matrix.  Each  of  the  matrix 
components  are  in  general  complex.  Since  the  wavefronts 
can  be  viewed  from  reference  points  at  the  immediate  input 
and  output  of  the  system,  2=0  and  z'  =0.  The  Jones  ma¬ 
trix  thus  presents  information  about  the  effect  of  a  par¬ 
ticular  optical  system  on  the  azimuthal  orientation  rela¬ 
tionship  between  the  input  and  output  reference  axes.  Any 
phase  modification  of  the  electric  field  components  caused 
by  the  system  are  contained  in  the  Jones  matrix. 

Another  matrix  of  primary  importance  in  optical  sys¬ 
tems  is  the  rotation  matrix.  This  matrix  allows  a  coordi¬ 
nate  rotation  of  the  Jones  vector's  reference  axes.  Figure 
2.  presents  the  graphical  aspects  of  this  rotation.  Shown 
here  is  the  rotation  of  an  electric  field  vector  with  ref¬ 
erence  axes  x,y  into  axes  labelled  x',y'.  The  angle  of 
rotation  is  7.  Using  the  Jones  vector  it  is  apparent  that 
the  rotation  can  be  described  in  the  following  manner: 

r£)('1_rc«4T  smTir^x"! 

UvJ  CosTfJie-jJ  (6) 

where  the  matrix  performing  the  rotation  operation  on  the 
Jones  vector  is  the  rotation  matrix.  A  counter-rotation 
from  x',y'  to  x,y  is  also  available  by  merely  changing  the 
sign  of  7. 

Armed  with  the  Jones  formalism,  the  theory  concerning 


the  measurement  of  the  polarization  ellipse  is  straight¬ 
forward  (2) .  The  PSA  system  and  its  reference  axes  are  il¬ 
lustrated  in  Fig.  3.  The  x  and  y  axes  are  taken  as  respec¬ 
tively  the  horizontal  and  vertical  axes  viewed  with  the 
beam  traveling  toward  the  observer.  The  electric  field  at 
the  PMT  input,  can  be  expressed  as  the  product  of  the 
following  matrices  operating  on  the  electric  field  exiting 
the  polarizer,  Ep^, 

Eo  =  "f;R(A^5R(-P)Ep, 

The  R()  matrices  represent  the  coordinate  rotations  refer¬ 
enced  to  the  polarizer  (P)  and  analyzer  (A)  angles. 

The  polarizer  allows  passage  of  light  linearly  polar¬ 
ized  in  the  x'  direction.  Therefore  the  effect  of  the 
polarizer  section  is  given  by  the  following 
counter-rotation : 


(8) 


The  sample  is  described  by  a  Jones  matrix  having  the 

complex  reflection  coefficients  r  and  r  as  matrix  el- 

p  s 

ements.  These  reflection  coefficients  contain  the  ratio  of 
the  reflected  to  incident  electric  field  and  are  given  by 


(9a) 

Here  the  subscripts  p  and  s 


f  _  Esr  ,  Elf 

and  ’S~  ■=— .  (9b) 

(from  the  German  "senkrecht”) 


refer  to  the  field  components  oriented  parallel  and  per¬ 
pendicular  to  the  plane  of  Incidence  of  the  sample.  The 


Sample. 

Suf  fa  CC. 

Figure  3.  The  PSA  system  showing  the  reference  axes  for 
each  optical  component. 
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subscripts  r  and  i  stand  for  reflected  and  incident.  The 
sample  alignment  procedure  prior  to  measurement  orients 
the  plane  of  incidence  so  that  the  p  and  s  axes  coincide 
with  the  X  and  y  reference  axes.  The  transfer  matrix  of 
the  sample  is  thus  given  by 


(10) 


The  analyzer  will  permit  the  transmission  of  the 
fleeted  electric  field  components  oriented  only  in  the 
direction.  Thus  a  rotation  from  the  x,y  to  x''  axes 
needed: 


TiRw) 


CflsA 

Sin  A 


Sin  A  "I 
Cos  A  i 


(11) 


re¬ 
x'  ' 


is 


The  PSA  system  transfer  function  is  obtained  by  com 
bining  equations  8,  10  and  11: 

CsnoirCssA  SiftA  "If^P  o'!  PcosP  -S'mP 
^  Lo*Jb^^A  cojaJLo  ’fjJLS'ftP  c«sP 

=  CosP  Cos  A  +  Vj  Sin?  Sift  A  ■] 

L  o  J  (12) 


The  intensity  at  the  PMT  input  is  represented  by  the  fol¬ 
lowing  proportionality; 

A  proportionality  is  used  since  the  characteristic  imped¬ 
ance  is  omitted.  The  expression  for  intensity  is  taken  out 
of  matrix  form  at  this  point  since  the  value  seen  at 
the  PMT  input  depends  only  upon  the  x''  orientation  of  the 


rotating  analyzer  fast  axis.  The  intensity  at  the  detector 
is  easily  determined  by  inserting  Ej^  and  from  equation 
12  into  equation  13  and  performing  the  mathematical  op¬ 
erations: 

Eq  i<.[<pC>iPcesA 

=  Ce«?A  ♦  Sin* A 

t(rprs’^+V»fs)Sih?  Ct.i?SiftAC6sA 

=  1  Ppl^Cos’P  (I  +C1.S  2A)  +  P  (l-Cos  2A) 

+  iRe  Sin P Cos  P  Sin  ZA 

=  |(V\^Cas*P  +l^s\*S!n*P+DtVl’Cos*P-lrs\'Si,?p]ices2A 

■^iReCPprs^iSinPCosP  SinZA 

Ip  can  be  normalized  so  that  its  has  an  average  value  of 
1,  thus  yielding 


Tp  oc  1  t  1  -  \ox\  P 

IHf-'TAn-P 

+  ^  C  )  "tftn  P  Sip(2A 

l^r+Tan-^P 


(14) 


This  expression  can  be  further  simplified  using  the  rela 
tionship  between  the  ellipsometric  parameters  V  and  A  and 
the  complex  reflection  coefficients  of  equation  9a  and  9b 


Ijj  then  becomes 

T  «1+ ^'r»n'/’CosA~TMLgs;^2A 

^  TftnV*T«»fP  lao’g- tToA^P 

(16) 

which  is  reducible  to 

oc  Cos  ZA  t  p  Sin  ZA  (17) 

Ip  (which  is  normalized  in  this  case)  is  a  sum  of 
sinusoidal  terms  riding  on  a  DC  level.  This  intensity  is 
the  quantity  measured  as  explained  in  the  section  on  the 
PSA  system.  The  coefficients  a  and  p  are  the  normalized 
Fourier  coefficients  of  the  intensity.  The  values  of  a  and 
P  are  determined  after  data  acquisition  by  a  Fourier 
analysis  of  Ip  contained  in  a  computer  program. 

To  summarize,  VASE  measurement  thory  is  based  on  the 
following  process: 

1)  Measure  the  light  intensity  coming  from  the  sample. 

2)  Determine  the  Fourier  coefficients  a  and  p. 

3)  Invert  the  Fourier  coefficients  to  find  Tan0  and  CosA 
(V>  and  A)  . 

4)  Determine  the  reflection  coefficients  r^  and  r  from  V 

p  s 

and  A. 

5)  The  dielectric  properties  ,  optical  constants,  alloy 


composition,  and  layer  thicknesses  are  calculated  using 
the  reflection  coefficients  and  the  angle  of  incidence. 
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Shown  in  Fig.  1  is  a  schematic  diagram  of  the  rotating 
analyzer  VASE  used  for  the  experimental  phase  of  this  re¬ 
search  effort.  The  system  is  based  upon  a  design  by  Aspnes 
and  Studna  (1).  However,  the  system  at  the  University  of 
Nebraska  has  been  modified  to  include  variable  angle  of 
incidence  capability.  The  system  is  basically  composed  of 
the  following  sections:  light  beam  generation, 
ellipsometer ,  and  data  acquisition  electronics. 

The  beam  generation  section  contains  the  lamp, 
monochromator,  and  collimator  optics.  The  lamp  is  a  75  W 
Xenon  short  arc  type  with  a  spectral  emission  in  the 
2,000  to  20,000  A  range.  The  arc  lamp  output  is  passed  to 
a  Kratos  GM252  monochromator  which  provides  a  narrow  band 
of  wavelengths  from  2500  to  8500  A  with  a  line  width  of 
approximately  20  A.  The  beam  divergence  at  the 
monochromator  output  is  controlled  by  an  off-axis 
parabolic  reflector  and  the  input  collimator  on  the 
polarizer. 

The  ellipsometer  is  a  Gaertner  model  L119  with  the 
capability  of  varying  the  angle  of  incidence,  <f>,  from  20* 
to  almost  90*  with  an  accuracy  of  +/-  0.01*  .  The 
monochromator  beam  enters  the  polarizer  and  exits  at  a 
known,  linear  polarization.  The  beam  strikes  the  sample 
becoming  elliptically  polarized  and  is  reflected  to  the 
rotating  analyzer.  The  rotating  analyzer  sweeps  the  ellip- 


tical  beam  at  3600  rpm.  Thus  a  signal  which  is  amplitude 
modulated  by  an  amount  determined  by  the  eccentricity  of 
the  polarization  ellipse  is  passed  to  the  photomultiplier 
tube  (PMT) . 

The  data  acquisition  section  converts  the  analog  PMT 
signal  voltage  to  a  digital  response  which  is  passed  to  a 
computer  for  analysis.  The  output  of  the  PMT  is  a 
sinusoidal  voltage  impressed  on  a  DC  level.  This  PMT  out¬ 
put  is  amplified  and  sent  to  the  A/D  converter.  The 
maximum  and  minimum  excursions  of  the  sinusoid  above  DC 
represent  respectively  the  long  and  short  axis  of  the  po¬ 
larization  ellipse.  A  reference  pulse  is  initiated  at  a 
calibrated  setting  of  the  analyzer.  The  first  count  pulse 
is  initiated  simultaneously  with  the  reference  pulse.  The 
signal  is  then  sampled  72  times  per  revolution  of  the 
analyzer.  The  sampling  period  is  variable  between  24  and 
300  cycles  per  reading.  Longer  sampling  periods  allow  in¬ 
creased  averaging  of  the  A/D  converter  output.  This  aver¬ 
aging  suppresses  noise  effects  in  spectral  regions  where 
the  signal  to  noise  ratio  is  low.  The  sampled  and  averaged 
signal  is  sent  in  digital  form  to  the  computer  where  a 
complex  analysis  yields  the  ellipsometric  parameters,  di¬ 
electric  (optical)  constants  as  a  function  of  wavelength. 
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The  experimental  data  taken  from  VASE  measurements  is 
modeled  using  the  Multiple  Angle  and  Wavelength  for  Effec¬ 
tive  Media  Approximation  (MAWEMA)  program  written  by  Mar¬ 
tin  C.  Rost  at  the  University  of  Nebraska.  The  overall 
purpose  of  the  program  is  to  fit  calculated  values  of  0 
and  A  to  the  experimental  values  of  the  two  ellipsometric 
parameters  as  a  function  of  wavelength.  MAWEMA  is  also 
capable  of  generating  rp  and  A  values  given  a  layered 
structure  with  fixed  layer  thicknesses  and  compositions. 

The  data  fitting  is  done  by  first  specifying  the  ap¬ 
propriate  variable(s)  (i.e.  layer  thickness,  alloy  compo¬ 
sition,  etc.)  in  the  model  file.  Using  initial  estimates 
of  those  variables,  V*  and  A  values  are  generated  using 
the  Fresnel  reflection  coefficients  for  a  multilayered 
structure.  The  calculated  tp  and  A  values  are  then  com¬ 
pared  with  the  experimental  ip  and  A  values  and  the  mean 
square  error  (MSE)  between  the  calculated  and  experimental 
values  is  determined.  The  standard  MSE  is  given  by 

MS  e  =  ^  I [CVi;  -  •I'.i  (A.;  - 

where  the  subscripts  c  and  i  stand  for  calculated  and  ex¬ 
perimental  respectively,  and  i  is  an  index  that  ranges 
over  all  the  m  measured  values  of  ip  and  A.  The  MSE  may 
also  be  determined  with  respect  to:  1)  V"  only,  2)  A  only, 
3)  TanV>  and  CosA,  4)  TanV»  only,  and  5)  Cos  A  only.  After 
MSE  determination,  the  constituent  variables  are  per- 


turbed  and  the  process  is  repeated.  Each  perturbation  is 
directed  toward  producing  the  smallest  MSE  possible 
between  the  calculated  and  experimental  values  of  the 
ellipsometric  parameters. 

Physically,  the  computations  are  carried  out  using 
three  files:  1)  the  data  input  file,  labeled  FILE. INN,  2) 
the  model  file,  labeled  FILE. MOD,  and  3)  the  data  output 
file,  labeled  FILE. OUT.  An  example  of  a  portion  of  an  INN 
file  obtained  from  a  GaAs/AlGaAs  structure  is  shown  in  Fig 
1.  This  file  is  produced  from  measured  data  that  is  trans¬ 
ferred  via  modem  from  the  VASE  laboratory  to  the  VAX  com¬ 
puter  system.  The  input  file  is  designated  as  either  EX¬ 
PERIMENTAL  or  GENERATED  on  the  first  line.  Only  the 
results  of  an  experimental  input  file  will  be  considered 
here.  Some  generated  results  are  provided  in  the  VASE 
studies  section  of  this  thesis  (M.C.  Rost's  manual  pro¬ 
vides  a  number  of  excellent  examples) .  The  first  column 
contains  the  wavelengths  of  the  measured  spectrum  (mea¬ 
sured  at  40  A  intervals  in  this  case) .  The  second  column 
contains  the  angle  of  incidence  with  the  third  and  fourth 
columns  manifesting  rj)  and  A  respectively.  The  MAWEMA  pro¬ 
gram  can  analyze  input  data  sets  with  single  or  multiple 
angles  of  Incidence. 

The  MOD  files  are  set  up  with  the  format  shown  in  Fig. 
2a.  An  example  model  file  for  a  GaAs/AlGaAs  multilayered 
sample  is  illustrated  in  Fig.  2b. 
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The  statistical  theory  behind  the  MAWEMA  program  is 
too  involved  and  voluminous  to  describe  here.  Therefore, 
the  explanation  of  the  MOD  structure  file  will  be  detailed 
only  with  respect  to  those  modeling  quantities  that  are 
alluded  to  in  the  latter  portions  of  this  thesis.  The  ex¬ 
planation  for  the  designations  in  the  MOD  file  format  of 
Fig.  2a.  is  as  follows: 

1.  EMA  Type  If  a  0  is  placed  in  this  position  the 
effective  medium  approximation  that  is  used  is 
Maxwel 1 -Garnett ;  a  1  signifies  Bruggeman  EMA.  The 
Maxwell -Garnett  and  Bruggeman  EMA  methods  are  used  to 
characterize  the  microstructure  of  film  layers  that  may 
consist  of  different  constituents  (i.e  a  roughness  layer 
that  is  a  mixture  of  50%  GaAs  and  50%  AlGaAs) .  The 
Maxwell-  Garnett  method  solves  for  dielectric  properties 
by  modeling  a  layer  as  a  host  material  containing  isolated 
"islands”  of  other  dielectric  materials.  The  Bruggeman 
method  assumes  a  random  mixture  of  a  number  of  dielectric 
materials.  The  Bruggeman  microstructure  is  thus  a  series 
of  interconnected  sections  of  differing  dielectric  materi¬ 
als.  Neither  of  these  EMAs  is  typically  used  in  the  solu¬ 
tion  of  the  layer  thicknesses  and  alloy  composition  of  MBE 
grown  samples  since  the  layer  interfaces  are  assumed 
abrupt  and  the  materials  in  the  layers  are  well  known.  In 
the  example  of  Fig  2b.  a  Bruggeman  EMA  is  designated  but 
is  not  required  in  the  model. 
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EXPERIMENTAL 
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Figure  1.  Example  of  an  INN  file  used  in  MAWEMA. 
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Figure  2b.  An  example  of  the  MOD  file  for  a  GaAs/AlGaAs 
structure. 


2.  MSE  TYPE  This  entry  is  an  integer  between  0  and  5 
which  determines  which  ellipsometric  parameter (s)  will  be 
minimized  in  the  analysis  as  explained  previously.  In  the 
example  the  MSE  is  being  calculated  with  respect  to  both  ^ 
and  A. 

3 .  ISEN  is  an  integer  entry  between  0  and  3  that 
specifies  whether  to  print  the  sensitivities  matrix  and/or 
the  90%  confidence  limits  for  each  variable.  The  90%  con¬ 
fidence  limits  are  based  on  Gaussian  statistics  and  yield 
a  measure  of  the  uncertainty  of  the  calculated  data  with 
respect  to  the  experimental  data.  Since  a  variable  value 
is  calculated  for  each  pair  taken  at  a  single  wave¬ 
length,  the  measured  semiconductor  system  is 
overdetermined.  The  final  calculated  values  of  layer 
thickness  and  composition  are  the  average  values  of  the 
variable  obtained  from  each  wavelength.  The  90%  confidence 
limits  are  the  variable  values  at  the  boundaries  of  90%  of 
the  area  underneath  a  Gaussian  bell  curve  centered  at  the 
average  value  of  that  variable.  In  the  example,  only  the 
90%  confidence  limits  will  be  placed  in  the  output  file. 

4.  ICOV  is  an  integer  entry  (0  for  "no",  and  l  for  ”yes") 
specifying  whether  or  not  to  print  the  covariance  matrix. 
From  the  covariance  matrix  the  correlation  matrix  is  cre¬ 
ated.  Correlation  is  the  degree  to  which  two  variables  are 
interdependent.  An  accurate  determination  of  correlation 
is  difficult  since  two  variables  may  have  a  low  correla- 
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tion  but  both  may  be  highly  correlated  to  a  third  vari¬ 
able.  High  correlation  exists  when  one  variable  is  changed 
in  the  model  and  the  calculated  values  for  that  variable 
and  another  variable  change  considerably.  Therefore 
cahnges  in  the  ellipsometric  data  cannot  be  directly  at¬ 
tributable  to  one  variable  only. 

5.  IPNT  is  an  integer  specifying  the  number  of  passes 
through  the  minimization  routine  for  which  the  intermedi¬ 
ate  results  are  to  be  printed.  If  the  MOD  file  of  Fig.  2b. 
is  used  the  results  of  every  20th  iteration  will  be  passed 
to  the  output  file. 

6.  ITMAX  is  an  integer  which  shows  the  number  of  it¬ 
erations  the  minimization  routine  will  perform  before  ter¬ 
minating  if  a  solution  hasn^t  been  found.  This  prevents 
MAWEMA  from  trying  to  minimize  the  MSE  indefinitely. 
Therefore,  99  or  less  iterations  will  be  performed  in  the 
example. 

7.  MSE VAR  is  a  real  value  which  indicates  by  what  percent 
a  variable  is  to  be  changed  in  order  to  determine  the 
change  in  MSE  from  one  iteration  to  the  next.  The  vari¬ 
ables  in  the  example  are  changed  by  1%  after  each  it¬ 
eration. 

8.  EPSl  is  a  real  number  which  specifies  at  what  relative 
value  of  MSE  from  one  iteration  to  the  next  the  program 
will  terminate.  MAWEMA  stops  when  the  difference  in  MSE 
between  subsequent  iterations  is  less  than  this  value.  The 
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minimization  will  stop  in  the  example  when  the  MSE  changes 
-4 

by  less  than  10 

9.  EPS 2  is  a  real  number  stating  the  amount  by  which  a 

variable  may  change  from  one  iteration  to  the  next  before 

the  model  fitting  process  is  ceased.  A  change  in  a  model 

-4 

variable  of  less  than  10  in  the  example  will  cause 
MAWEMA  to  terminate. 

10.  FILM  NUMBER  is  an  integer  value  which  labels  the  film 
layers.  The  lowest  numbered  film,  0,  is  the  bottom  layer 
in  the  sample.  In  the  example  there  are  four  films.  The 
bottom  film  is  a  GaAs  substrate  followed  by  an  AlGaAs 
layer,  a  GaAs  layer,  and  a  surface  layer  of  GaAs  native 
oxide. 

11.  DATA  FILE  NAME  calls  a  computer  data  file  which  con¬ 

tains  either  the  components  of  the  complex  permittivity 
(«!,  <2)  optical  constants  (n,  k)  over  the  wave¬ 

length  range  of  interest  for  that  particular  layer.  MAWEMA 
uses  the  optical  constants  in  calculating  the  Fresnel  re¬ 
flection  coefficients.  In  the  example,  the  optical  con¬ 
stants  for  the  following  materials  are  being  used:  GaAs 
(GAAS.DAT),  AIq  ^Ga^  gAs  (ALGAAS5.DAT),  Al^  ^Ga^  ^As 
(ALGAAS4.DAT),  Al^  7^®  (ALGAAS3.DAT),  AIq  2^30  8^® 

(ALGAAS2.DAT),  Al.  .Ga.  _As  (ALGAAS1.DAT).  The  optical 

U  •  X  U  •  7 

constants  in  the  .DAT  files  have  been  obtained  from  pub¬ 
lished  data  (1,2,3)  (some  measured  using  ellipsometry) . 
The  AlGaAs  materials  with  different  composition  assigned 


to  layer  1  are  used  in  a  subroutine  that  provides  interpo¬ 
lated  values  of  n  and  k  taken  from  the  .DAT  files.  This 
allows  the  composition  value  to  be  varied  between  x  =  0 
and  X  =  1. 

12.  FRACTION  states  what  decimal  percentage  of  that  layer 
is  the  constituent  in  the  data  file.  This  percentage  is 
used  in  the  Maxwell -Garnett  and  Bruggeman  EMAs.  The  0.3  in 
the  example  means  that  layer  1  is  AlGaAs  with  x  =  0.3. 

13 .  IFRAC  is  an  integer  stating  whether  the  FRACTION  is 
variable  (0) ,  fixed  (1) ,  or  dependent  upon  the  variable 
value  (2) .  The  AlGaAs  composition  in  layer  1  is  then  vari¬ 
able. 

14.  THICKNESS  (real  variable)  is  the  thickness  of  a  par¬ 
ticular  layer  in  angstroms.  In  the  GaAs/AlGaAs  .MOD  file 
shown,  the  bottom  layer  is  GaAs  substrate  with  a  thickness 
that  is  optically  infinite,  therefore  no  layer  thickness 
is  provided.  The  AlGaAs  layer  is  800  A  thick,  the  GaAs 
layer  500  A  and  the  oxide  layer  is  20  A. 

15.  ITHICK  is  an  integer  indicating  whether  the  layer 
thickness  is  variable  (0)  or  fixed  (1) .  All  layer  thick¬ 
nesses  in  the  example  are  variable  except  the  substrate. 

16.  I PEP  is  the  inter-film  dependence  integer.  Each  con¬ 
stituent  is  assigned  a  dependence  number  indicating 
whether  the  film  has  the  same  material  in  it  as  another 
film.  The  Interdependent  films  in  the  example  .MOD  file 
are  the  GaAs  films,  meaning  that  the  same  optical  con- 
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stants  are  used  for  each  of  the  GaAs  layers. 

17.  lose  is  an  integer  designating  whether  the  optical 
constants  from  a  .DAT  file  or  an  oscillator  data  file  is 
to  be  used.  An  oscillator  data  file  is  used  if  the  optical 
constants  are  expressable  as  Lorentzian  oscillators.  Each 
oscillator  has  a  resonant  frequency  (wavelength) .  Oscilla¬ 
tors  can  be  used  if  peaks  are  observed  in  the  experimen¬ 
tal,  optical  constant  data.  Single  or  multiple  oscillators 
can  then  be  used  to  fit  the  data  by  placing  an  oscillator 
at  the  wavelength  of  the  peak(s) .  No  oscillators  are  used 
in  the  example. 

18.  n _ k  These  are  the  number  of  n-k  data  sets  to  be 

solved  for  using  MAWEMA.  The  integer  0  means  that  they  are 
variable  and  that  n  and  k  are  to  be  determined  using  the 
experimental  dat  .  The  integer  1  is  used  for  fixed  values. 
The  number  of  n-k  sets  must  correspond  to  the  number  of  0 
and  A  experimental  data  sets.  In  this  case  the  sample  was 
measured  at  20  different  wavelengths  and  the  fixed  optical 
constants  at  each  of  those  wavelengths  must  be  called  up 
from  the  .DAT  files.  Notice  that  the  n-k  sets  do  not  need 
to  be  repeated  for  films  that  are  interdependent  since  the 
optical  constants  for  that  film  have  already  been  used. 

It  should  be  stressed  that  the  statistical  parameters 
determined  by  the  MAWEMA  program  (i.e.,  MSE,  90%  confi¬ 
dence  limits,  and  correlation  )  are  not  singularly  in¬ 
dicative  of  the  correctness  of  the  data  fitting  procedure. 
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This  is  due  to  the  fact  that  the  results  are  dependent 
upon  the  initial  model  used.  When  modeling  MBE  grown  GaAs 
samples,  the  values  of  layer  thickness  and  composition 
provided  by  the  growers  are  considered  to  be  at  least 
"ballpark”  figures.  The  statistical  values  are  implemented 
mainly  for  the  comparison  of  different  minimization  runs. 
A  good  modeling  procedure  is  to  use  different  starting 
values  of  layer  thickness  and  composition  and  compare  the 
results  of  each  analysis.  If  the  layer  thicknesses  and 
compositions  are  the  same,  then  an  accurate  analysis  is 
assured. 

The  output  file  (.OUT)  contains  the  following: 

1)  The  experimental  and  calculated  values  of  rji  and  A  at 
the  wavelengths  measured. 

2)  The  experimental  and  calculated  dielectric  constants  of 
each  layer.  The  dielectric  constants  are  related  to  the 
ellipsometric  parameters  by 

(2) 

The  dielectric  constants  can  be  expressed  in  terms  of  the 
optical  constants  by 

»  n’--  6,  =  Zr>A 

*  (3) 

The  index  of  refraction  of  an  absorbing  material  is  a  com- 
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plex  quantity  given  by 

Ni  =T»  +  lA 

(4) 

This  value,  N,  is  used  in  the  Fresnel  coefficients  and 
provides  the  phase  difference  between  the  reflection  coef¬ 
ficients  and  hence  elliptical  polarization  of  the  re¬ 
flected  beam. 

3)  Minimization  statistics,  i.e.,  the  MSB,  90%  confidence 
limits  and  the  correlation  matrix. 
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Chapter  2. 

The  Physics  and  Technology  of  GaAs/AlGaAs  Heterostructures 


A.  MODFET  Technology 

The  advent  of  the  layered  semiconductor  super lattice, 
first  proposed  by  Esaki  and  Tsu  in  the  late  1960s,  led  to 
the  development  of  the  first  modulation-doped  field  effect 
transistor  (MODFET)  by  Dingle,  Stormer  and  Gossard  at  Bell 
Laboratories  in  1978.  Since  then,  improvements  in  device 
speed,  power  consumption  and  operating  frequencies  have 
occurred  through  the  further  development  of  MODFET  tech¬ 
nology. 

MODFETs  are  formed  by  growing  alternate  layers  of 
doped  AlGaAs  on  undoped  GaAs  (hence  the  term 
"modulation-doped”) .  An  energy  band  diagram  for  the 
MODFET  is  shown  in  Fig.  1.  Since  the  GaAs  has  a  smaller 
band  gap  energy  than  the  Alj^Gaj^__j^As,  a  band  discontinuity 
(dependent  on  the  AlGaAs  mole  fraction  x)  exists  at  their 
juncture  and  quantized  electron  and  hole  levels  occur  in 
the  GaAs  layer.  The  donor  doping  level  in  the  AlGaAs  is 
at  a  higher  energy  than  the  conduction  band  energy  of  the 
GaAs.  Consequently,  the  electrons  in  the  AlGaAs  diffuse 
into  the  GaAs  layer  until  charge  equlibrium  occurs.  The 
carrier  diffusion  produces  a  large  electric  field  (-  100 
kV/cm)  which  causes  band  bending.  The  band  bending  re¬ 
sults  in  the  formation  of  a  quasi-triangular  potential 
well.  The  quantization  takes  place  in  the  dimension  per¬ 
pendicular  to  the  heterostructure,  resulting  in 


unrestrained  electron  motion  parallel  to  the  interface 


(1) .  The  electron  system  is  thus  two  dimensional  and  is 
referred  to  as  a  two-dimensional  electron  gas  (2DEG) .  The 
separation  of  ionized  impurities  in  the  AlGaAs  from  the 
conduction  electrons  in  the  GaAs  creates  increased  carrier 
mobility  due  to  the  reduction  of  ionized  impurity  scatter¬ 
ing. 

A  typical  MODFET  device  is  illustrated  in  Fig.  2  (3). 
This  device  is  employed  in  a  normally-off  configuration 
where  the  gate  depletion  region  extends  into  the 
two-dimensional  electron  gas  connecting  the  source  and 
drain,  cutting  off  conduction.  In  the  normally  off  condi¬ 
tion,  the  gate  and  source  are  kept  at  the  same  voltage. 
Once  the  gate  voltage  is  0.6  to  0.8  V  higher  than  the 
source,  the  depletion  region  is  moved  up  into  the  undoped 
AlGaAs  layer  and  conduction  ensues.  The  isolation  region 
serves  to  separate  MODFET  devices  on  the  same  integrated 
circuit.  The  n^  GaAs  cap  layer  prevents  oxidation  of  the 
AlGaAs  layer  and  allows  low  resistance  ohmic  contacts  for 
the  drain  and  source  (2) .  The  ohmic  metal  is  usually  n 
AuGe/Ni/Au  alloy.  A  Schottky  contact  is  formed  when  the 
gate  metal  is  placed  on  top  of  the  AlGaAs  layer. 
Typically,  a  thick,  undoped  GaAs  buffer  layer  (~  1  urn) 
is  placed  between  the  MODFET  and  the  substrate  to  provide 
a  smooth  surface  on  which  to  grow  the  MODFET.  The  sub¬ 
strate  is  composed  of  GaAs  which  is  Cr-doped.  This  pro¬ 
vides  a  semi-insulating  layer  which  prevents  stray  conduc- 


Figure  1.  MODFET  energy  band  structure  showing  band  offset 
AE^  and  AE^,  and  quantized  energy  levels  Eq  and  E^^  in  the 
quasi-triangular  potential  well.  (From  Ref.  2) 


Figure  2.  Schematic  diagram  of  MODFET  construction  for  ic 
application.  (From  Ref.  3) 


tlon  and  decreases  parasitic  capacitance. 

From  a  device  standpoint,  the  main  advantage  of  the 

MODFET  is  the  positioning  of  a  large  electron  density 
19  —3 

(-10  cm  )  in  a  thin  layer  (-80  A)  very  close  to  the 
gate  (-300  A) .  These  factors  permit  large  currents  (-500 
mA/mm)  to  be  switched  at  very  high  speeds  (-lo”^®  s  @  300 

K)  allowing  exceptionally  fast  capacitor  charging  times 
with  reduced  power  consumption  (-1  mW  @  300  K)  (4) . 

To  enhance  the  device  properties  of  the  MODFET, 
undoped  AlGaAs  spacer  layers  and  selective  doping  are 
used.  The  spacer  layer  is  placed  between  the  doped  AlGaAs 
and  undoped  GaAs  layers.  The  spacer  layer  (shown  sche¬ 
matically  in  Fig.  3)  provides  increased  isolation  between 
the  ionized  donors  in  the  doped  AlGaAs  and  the  2 DEG.  This 
quality  further  reduces  ionized  impurity  scattering  and 

enhances  the  band  structure  by  increasing  the  conduction 

.  .  12 

band  offset.  An  planar  layer  of  single  atom  Si  (-  10 

—2 

cm  )  -  called  "delta”  doping  -  can  be  placed  in  the 

undoped  spacer  layer  close  to  the  GaAs  well.  The  close 
proximity  of  the  delta  doped  layer  to  the  well  produces  an 
additional  supply  of  conduction  electrons,  thus  increasing 
the  2DEG  concentration.  Improvements  in  2DEG  concentra¬ 
tion  of  30%  with  almost  the  same  mobility  have  been  real¬ 
ized  using  delta  doping  (5) . 

MODFET  structures  where  the  GaAs  layer  is  grown  on 
top  of  the  AlGaAs  have  received  much  attention  in  recent 


Ns  =  7  X  10"  cm"^  {300K) 


Fermi  Level 


Figure  3 .  MODFET  energy  band  diagram  showing  the  undoped 
AlGaAs  spacer  layer  and  depletion  region.  (From  Ref.  6) 


n*  GaAs 


AlasGaojAs 
Spacer - 1 


n”  GoAs 


n.AlojGoojAs 


Buffer 


GoAs  Buffer 


S.I.GoAs  Substrofe 


Figure  4.  Schematic  of  an  "inverted”  MODFET  (from  Ref.  4) 
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years.  A  typical  "inverted”  configuration  is  provided  in 
Fig.  4.  It  is  composed  of  an  undoped  AlGaAs  buffer 
layer,  a  doped  AlGaAs  layer,  an  undoped  AlGaAs  spacer 
layer  and  a  GaAs  cap  layer.  The  GaAs  layers  at  the  source 
and  drain  are  doped  n^  in  order  to  make  ohmic  contacts  at 
these  points. 

One  problem  with  this  technology  is  that  inverted 
structures  have  exhibited  degraded  mobility  character¬ 
istics.  A  number  of  growth  techniques  have  been  employed 
to  counter  this  undesirable  quality  and  are  discussed 
later. 

An  inverted  MODFET  using  a  wide  bandgap  superlattice 
is  shown  in  Fig.  5.  The  25  period  undoped  superlattice  on 
top  of  the  GaAs  buffer  layer  prevents  the  transfer  of 
electrons  into  the  buffer  layer.  The  10  period 
superlattice  contains  a  15  A  donor  doped  (10  cm  )  GaAs 
layer  sandwiched  between  undoped  5  A  GaAs  layers.  The 
dopants  reside  at  shallow  donor  levels  in  the  GaAs  portion 
of  the  superlattice  and  provide  higher  electron  sheet 
densities  for  the  2DEG.  Thus  this  superlattice  takes  the 
place  of  the  doped  AlGaAs  layer  found  in  the  conventional 
MODFET.  The  5  A  undoped  GaAs  layers  in  this  superlattice 
are  used  to  prevent  Si  incorporation  in  the  AlAs.  The  four 
period  undoped  superlattice  provides  a  spacer  region  which 
further  reduces  ionized  impurity  scattering  between  donors 
in  the  doped  superlattice  and  the  2DEG. 
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Figure  5.  Cross  sectional  view  of  an  inverted  MODFET  with 


GaAs/AlAs  superlattices  (from  Ref.  4) . 
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Figure  6.  A  double-heterointerface  MODFET  grown  on  a  three 


superlattice  structure  (from  Ref.  4) . 
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A  third  type  of  MODFET  employs  quantum  wells  and  uses 
both  the  inverted  and  a  normal  interface.  A  thin  GaAs 
layer  (~  300  A)  is  sandwiched  between  two  doped  AlGaAs 
layers  forming  a  quantum  well.  The  2DEG  concentration  is 
increased  since  electrons  are  available  from  both  barrier 
layers  at  each  interface,  thereby  increasing  the  current 
carrying  capability  (4) .  The  gate  to  2DEG  distance  is  kept 
small  since  the  GaAs  layer  is  thin.  An  example  of  a 
double-heterointerface  (quantum  well)  MODFET  is  given  in 
Fig.  6.  The  structure  has  a  1.1  urn  GaAs  buffer  layer  with 
a  33  period  undoped  GaAs  superlattice,  a  10  period  n-doped 
GaAs  superlattice,  and  a  one  and  a  half  period 
superlattice  used  as  a  spacer  layer.  The  10  period 
superlattice  is  used  as  a  carrrier  supplying  layer  for  the 
quantum  well  similar  to  that  in  Fig.  5.  The  33  period 
superlattice  is  used  as  barrier  layer  to  prevent  electron 
transfer  into  the  buffer  layer. 

MODFET  integrated  circuits  for  high  speed  computer 
applications  are  under  development  at  numerous  laborato¬ 
ries.  One  example  is  an  integrated  frequency  divider  cir¬ 
cuit  built  with  0.5  urn  gate  length  MODFETs  (7).  Perfor¬ 
mance  at  5.5  GHz  at  300  K  has  been  demonstrated  with  logic 
delays  of  22  ps/gate  and  power  dissipation  of  2.9  mW/gate 
and  a  fanout  of  two.  The  speed  of  this  device  is  roughly 
three  times  faster  than  that  of  conventional  GaAs  MESFET 
technology.  Of  critical  importance  is  the  layer  thickness 
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of  the  n-AlGaAs  layer  below  the  gate.  This  thickness  de¬ 
termines  threshold  voltage  and  transconductance.  Wafer 
homogeneity  is  thus  crucial  and  must  be  carefully 
monitored.  VASE  has  been  implemented  to  fufill  the  task 
of  determining  the  thickness  uniformity  in  MBE  growth  (8) . 

As  another  example,  a  1-kbit  static  RAM  has  been  de¬ 
veloped  using  MODFET  LSI  technology  (9) .  This  IC  is  fully 
operational,  employing  7,244  MODFETs  on  a  chip  with  dimen¬ 
sions  3.0  mm  X  2.9  mm  and  100  chips  per  wafer.  Deviations 
in  access  time  were  +/“  4%  which  is  much  smaller  than  that 
of  GaAs  MESFET  SRAMs.  Uniformity  in  gate  threshold 
voltage  was  excellent  (with  variations  less  than  +/-  0.05 
V)  due  to  the  high  layer  thickness  uniformity  (+/”  1%)  ob¬ 
tained  during  MBE  growth  and  a  self-terminating  selective 
dry  etching-process. 
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B.  Molecular  Beam  Epitaxial  Growth  of  Quantum  Wells  and 

Superlattices 

Molecular  beam  epitaxy  (MBE)  is  a  method  for  growing 
very  thin  semiconductor  layers  of  high  quality  and  unifor¬ 
mity.  Due  to  the  lower  growth  temperatures  used  in  the 
MBE  process  (-  625  *C) ,  more  abrupt  layer  interfaces  can 
be  grown  when  compared  with  the  vapor  phase  epitaxy  (VPE) 
(~  700 *C)  and  liquid  phase  epitaxy  (LPE)  (  -  800 *C)  meth¬ 
ods  (1) .  Metal  organic  chemical  vapor  deposition  (MOCVD) 
is  a  technology  with  growth  temperatures  and  defect  densi¬ 
ties  comparable  to  MBE.  MOCVD  has  a  higher  wafer  through¬ 
put  and  lower  basic  system  cost.  Considerable  research 
and  development  effort  is  going  into  the  improvement  of 
MOCVD.  However,  the  high  degree  of  growth  control  obtain¬ 
able  with  the  MBE  shutters  allows  the  MBE  system  to  pro¬ 
duce  the  thinnest  layers  and  most  abrupt  dopant  transi¬ 
tions  (2)  at  the  present  time.  The  above  mentioned 
advantages,  combined  with  the  capability  of  monitoring 
the  wafer  surface  during  growth,  make  MBE  an  ideal  instru¬ 
ment  for  realizing  the  structural  properties  required  for 
quantum  wells  and  superlattices. 

The  MBE  process  is  one  of  species  evaporation  under 
ultra-high  vacuum  (UHV) ,  involving  the  reaction  of  thermal 
atomic  (or  molecular)  beams  with  a  heated  substrate.  The 
growth  rate  is  generally  low  -  around  1  mm/hr  (i.e.  ap¬ 
proximately  one  monolayer  per  second) .  This  rate  allows 
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Figure  2.  MBE  growth  system  layout  (from  Ref.  1) 


cient  of  the  AS2  is  unity.  Theoretically,  each  Ga  atom 
site  will  have  an  AS2  dimer  attached  to  it  and  complete 
coverage  of  the  Ga  monolayer  by  an  As  monolayer  will  take 
place.  If  the  AS2  coverage  is  not  complete,  a  Ga  rich  sur¬ 
face  is  produced  along  with  an  increased  defect  density 
producing  deep  impurity  levels. 

A  typical  MBE  system  is  shown  in  Fig.  2.  The  wafer  is 
placed  in  the  MBE  machine  through  the  cassette  entry  lock. 
From  there,  it  is  transported  to  the  preparation/analysis 
chamber  where:  1)  surface  analysis  using  Auger  electron 
spectroscopy  (AES)  or  X-ray  photoelectron  spectroscopy 
(XPS)  is  performed,  2)  the  sample  is  preheated  in  order  to 
outgas  it,  and  3)  metallic  contacts  are  deposited.  The 
sample  is  then  sent  to  the  MBE  chamber  where  the  molecular 
beams  are  generated  from  Knudsen  effusion  cells  and  the 
fluxes  are  measured  by  a  mass  spectrometer.  Several 
Knudsen  cells  may  be  incorporated  into  the  growth  chamber 
in  order  to  dope  the  wafer  or  provide  multicomponent  al¬ 
loys  and  compounds. 

The  growth  process  can  be  monitored  by  surface  stud¬ 
ies  using  reflection  high  energy  electron  diffraction 
(RHEED) .  The  RHEED  system  measures  surface  quality  by  di¬ 
recting  an  electron  beam  on  the  sample  at  a  glancing 
angle.  The  reflected  beam  produces  diffraction  patterns  on 
a  fluorescent  screen.  The  relative  diffraction  intensi¬ 
ties  indicate  monolayer  growth  with  the  maximum  intensity 


occurring  for  each  successive  monolayer.  Due  to  the  lack 
of  total  coverage  of  successively  grown  layers,  MBE  sur¬ 
faces  become  less  smooth  as  growth  continues.  The  contri¬ 
butions  to  the  diffraction  pattern  from  these  partiall'' 
grown  layers  (unsmooth  layers  produce  diffraction  minima) 
causes  the  RHEED  intensities  to  die  out  after  ap¬ 
proximately  100  monolayers.  This  limitation  renders  the 
RHEED  process  an  incomplete  method  for  monitoring  layer 
growth,  and  provides  an  ideal  application  for  VASE  which 
can  be  applied  in  situ  or  after  growth. 

Of  particular  interest  in  MBE  growth  is  impurity  in¬ 
corporation  during  the  growth  process.  When  the  binary 
(GaAs)  layer  is  grown  on  top  of  the  ternary  (AlGaAs) ,  the 
interface  has  been  observed  to  contain  more  impurities  and 
roughness  than  with  the  ternary  grown  on  binary  (6,7,8). 
This  effect  is  prevalent  in  superlattices  which  contain 
multiple  interfaces  and  "inverted”  MODFETS.  Decreased 
mobility  is  observed  (compared  to  the  normal  interface) 
due  to  the  presence  of  acceptor  levels  in  the  GaAs  pro¬ 
duced  by  growth-related  carbon  impurities  (7) .  The  free 
electrons  recombine  with  the  acceptor  level  holes,  thus 
reducing  the  number  of  electrons  available  for  conduction. 

Roughness  at  the  interface  has  been  attributed  to:  1) 
noncrystallinity  associated  with  slight  lattice  mismatch, 
2)  differences  in  migration  length  (the  distance  an  atomic 
or  molecular  species  travels  along  the  surface  before  in- 


corporation)  between  Ga  and  Al,  and  3)  residual  oxygen  in¬ 
corporation.  The  lattice  mismatch  between  GaAs  and  AlGaAs 
(~  0.1  %)  can  lead  to  tensile  strain  perpendicular  to  the 
interface  causing  roughness  (9) .  Roughness  on  the  order 
of  multiple  monolayers  is  possible  due  to  the  larger  mi¬ 
gration  length  of  Ga  (150  A)  than  A1  (35  A)  from  the 
AlGaAs  surface  during  growth.  The  surface  migration  rate 
for  Ga  atoms  is  longer  than  that  of  the  A1  atoms  due  to 
the  Ga-As  bonds  being  weaker  than  the  Al-As  bonds.  Conse¬ 
quently,  the  inverted  interface  will  be  rougher  than  the 
normal  interface  where  only  GaAs  atoms  are  present  (10) . 
The  high  reactivity  of  oxygen  with  Al  is  also  a  possible 
cause  of  impurity  related  roughness.  Surface  strain  on  the 
growing  AlGaAs  layer  induces  oxygen  precipitation.  Subse¬ 
quently,  oxygen  burial  takes  place  in  the  first  few 
monolayers  of  the  GaAs  with  the  O  atoms  being  substitu- 
tionally  incorporated  at  As  sites  (8) .  Also,  oxygen  in¬ 
corporation  contributes  to  deep  levels  within  the  GaAs, 
causing  undesirable  recombination. 

Several  growth  techniques  have  been  implemented  to 
counter  the  problems  associated  with  inverted  interfaces. 
Residual  impurities  (such  as  O  and  C)  float  on  the  AlGaAs 
growth  surface  and  prevent  lateral  propagation  of  the 
source  species  thus  causing  nonuniform  surface  growth  and 
roughness.  The  growth  of  a  superlattice  below  the  GaAs 
quantum  well  has  reportedly  yielded  considerably  smoother 


GaAs-AlGaAs  interfaces  (6,9).  This  improvement  is  attrib¬ 
uted  to  lattice-strain  induced  diffusion  of  carbon  and 
oxygen  impurities  into  the  AlGaAs  layers  of  the 
superlattice  during  cooling  to  room  temperature.  A  sec¬ 
ond  technique  involves  the  interruption  of  the  growth  af¬ 
ter  the  deposition  of  the  AlGaAs  layer.  This  reduces  the 
roughness  to  only  a  single  monolayer  at  the  AlGaAs  surface 
and  reduces  residual  impurity  diffusion  in  the  first  few 
GaAs  monolayers  (10) . 

Improvements  in  MBE  methods  have  permitted  the  growth 
of  "strained”  layer  superlattices.  These  structures  are 
characterized  by  materials  with  lattice  mismatches  in  the 
0.1%  to  1%  range  which  greatly  enhances  the  number  of  ma¬ 
terials  that  can  be  grown  epitaxially  on  a  given  substrate 
(11) .  For  example,  InAlAs/InGaAs  MODFETs  (lattice  mis¬ 
match  approximately  2%)  grown  by  MBE  have  exhibited  room 
temperature  mobilities  50%  higher  than  similar  AlGaAs/GaAs 
systems  (12) .  Lattice  mismatch  between  the  layers  is  to¬ 
tally  accomodated  by  strain  in  the  layers  instead  of  at 
the  interfaces  where  strain  causes  dislocations  (13).  The 
electronic  properties  of  lattice  mismatched  layers  are 
different  from  those  of  bulk  materials  due  to  strain  ef- 
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Photon  absorption  in  a  semiconductor  can  create  a 
bound  electron-hole  (e-h)  pair  called  an  exciton.  The 
e-h  pair  is  bound  by  the  Coulombic  interaction  between 
the  two  particles  and  can  transport  energy  within  the 
crystal  without  transporting  charge  (1) .  Spectro¬ 
scopically,  these  excitations  are  exhibited  in  bulk  mate¬ 
rials  at  energies  equal  to  the  band  gap  of  the  semiconduc¬ 
tor  minus  the  binding  energy  of  the  exciton. 

Excitons  are  normally  classified  into  two  types: 
Frenkel  or  Mott-Wannier.  Excitons  belonging  to  the 
Frenkel  class  are  tightly  bound  e-h  pairs  localized  near  a 
host  atom.  They  appear  as  an  excited  state  of  a  single 
atom  which  may  move  from  site  to  site  on  the  lattice 
through  nearest  neighbor  coupling.  Frenkel  excitons  are 
prevalent  in  relatively  ionic  materials  and  hence  have 
large  binding  energies  (typically  1  eV) .  Mott-Wannier 
excitons,  on  the  other  hand,  are  weakly  bound  and  are 
present  in  covalently  bonded  materials.  The  Mott-Wannier 
exciton  is  thus  found  predominantly  in  semiconductors. 
The  Mott-Wannier  exciton  has  an  e-h  separation  distance 
which  is  large  compared  with  the  crystal  lattice  constant, 
producing  relatively  small  binding  energies  (typically 
from  3  to  40  meV)  (1,2). 

Excitonic  transitions  in  quantum  wells  and 
superlattices  are  easily  discernible  in  optical  spectra. 
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Since  electrons  and  holes  can  only  occupy  quantized  energy 
levels  within  a  quantum  well,  the  excitons  in  quantum 
wells  are  quantized.  The  quantized  recombination 
transitions  in  optical  spectra  produce  a  sharper  structure 
when  compared  to  transitions  from  the  valence  band  to  the 
conduction  band  continuum  in  a  bulk  semiconductor  (3). 
Consequently,  excitons  can  produce  changes  in  optical 
spectra  that  are  larger  than  those  for  equal  densities  of 
free  e-h  pairs  (i.e.,  the  exciton  resonance  peaks  are  dis¬ 
placed  from  the  continuum)  (4) . 

Excitonic  effects  are  enhanced  in  thin  layer  semi¬ 
conductors  due  to  increased  carrier  confinement.  The 
exciton  binding  energy  is  larger  in  a  quantum  well  than  in 
the  bulk.  A  decreased  physical  separation  between  the 
electron  and  hole,  as  a  result  of  carrier  confinement, 
causes  an  increase  in  the  Coulombic  potential  and  hence  a 
larger  binding  energy.  For  example,  =  4.2  meV  in  bulk 
GaAs  and  E^^^  =  9.1  meV  in  a  single,  50  A  width,  GaAs  quan¬ 
tum  well  with  an  Al  barrier  composition  x  =  0.3  (5).  The 
binding  energy  of  an  exciton  (either  in  a  bulk  material  or 
a  quantum  well  system)  is  described  by 

^=^5^  CTr)-'/2.;"  (1) 


where  p  is  the  effective  mass  of  the  exciton  and  n  is  the 


quantum  state. 

Of  particular  interest  are  electro-optical  effects 
produced  by  excitons  when  an  electric  field  is  applied  to 


a  quantum  well  or  superlattice  (6).  The  presence  of  the 
electric  field  promotes  a  high  degree  of  optical 
nonlinearity  in  these  structures  due  to  the  spectral 
shifting  of  the  absorption  edge  as  explained  below.  In 
addition,  excitonic  transitions  in  quantum  wells  take 
place  at  room  temperature  as  opposed  to  bulk  materials 
where  they  occur  only  at  very  low  temperatures.  This  dif¬ 
ference  between  bulk  and  quantum  well  excitons  occurs  be¬ 
cause  quantized  energy  levels  in  the  quantum  wells  pro¬ 
vide  much  sharper  excitonic  transition  energies  than  bulk 
material  even  at  room  temperature.  In  the  bulk  regime,  as 
the  temperature  is  increased,  electrons  are  thermal! zed  in 
greater  numbers  and  occupy  a  larger  number  of  states  in 
the  conduction  band  thus  producing  a  quasi-continuum  of 
exciton  energies  near  the  band  gap  energy.  Therefore,  bulk 
excitonic  transitions  are  thermally  broadened  at  room  tem¬ 
perature  and  do  not  possess  the  sharp  absorption  edges 
necessary  for  opto-electronic  applications.  This  combina¬ 
tion  of  nonlinearity  and  ambient  temperature  operation 
permits  the  application  of  these  structures  as 
electro-optical  switches  and  modulators. 

When  an  electric  field  is  applied  perpendicular  to  the 
layers,  a  small  broadening  and  shifting  of  the  exciton 
peaks  to  lower  energies  is  observed  -  commonly  referred  to 
as  the  quantum  Stark  effect  (6,7,8).  This  can  be  seen  in 
Fig.  1  where  the  peaks  in  the  absorption  spectrum  of  the 
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quantum  well  are  seen  to  move  to  lower  energies  with  the 
application  of  electric  field  strengths  of  16  kV/cm  in 
(b)  and  50  kV/cm  in  (c) .  The  shift  to  lower  energies  is 
caused  by  binding  energy  decreases  due  to  the  decrease  in 
Coulomb  attraction  between  the  e-h  pair  as  they  move  to 
opposite  sides  of  the  well.  In  Fig  1  (c) ,  the  excitonic 
peaks  are  much  less  resolvable  than  in  1  (a) .  This  broad¬ 
ening  of  the  peak  intensities  occurs  because  the  increase 
in  field  strength  promotes  tunneling  due  to  the  shifting 
of  the  carrier  concentration  and  a  decrease  in  the  barrier 
width.  This  process  is  shown  in  the  energy  band  diagram 
above  the  figure.  Therefore,  the  Stark  effect  produces 
nonlinear  absorption  characteristics  since  a  change  in  the 
electric  field  does  not  yield  the  same  absorption  spectrum 
as  the  no  field  condition.  This  property  is  an  advantage 
in  that  the  absorption  edge  can  be  moved  into  a  region 
where  the  sample  is  transparent  in  the  absence  of  a  field 
thus  making  it  possible  to  control  (modulate  or  switch) 
the  transmitivity  of  a  quantum  well  device. 

The  spectral  shift  is  also  a  function  of  well  width 
(8) .  In  quantum  wells  of  width  greater  than  50  A  the 
Stark  shift  is  prevalent  in  photoluminescence  (PL)  spec¬ 
tra.  However,  for  wells  less  than  50  A  wide,  the  exciton 
absorption  peaks  broaden  and  shift  to  higher  energies  when 
an  electric  field  is  applied.  For  wells  less  than  50  A 
the  Stark  effect  is  less  prevalent  because  the  electron  to 
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Figure  1.  Absorption  spectrum  for  a  multiple  quantum  well 
structure  with  applied  fields:  (a)  no  field, 
(b)  16  kv/cm  and  (c)  50  kV/cm  (from  Ref.  6) . 
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hole  distance  doesn't  increase  appreciably  with  bias. 
Therefore,  the  binding  energy  doesn't  noticeably  decrease 
with  an  applied  bias.  The  shift  to  higher  energies  of  the 
PL  peaks  of  the  thin  wells  is  suggested  as  being  caused  by 
an  anti-Stokes  electron-phonon  interaction.  These  elec¬ 
trons  thus  have  a  higher  recombination  energy.  The 
anti-Stokes  effect  occurs  for  the  thicker  wells;  however 
the  Stark  effect  is  stronger  for  these  wells.  These  two 
effects  are  shown  in  the  PL  spectra  of  Fig.  2b  for  a  mul¬ 
tiple  guantum  well  structure  (Fig  2a)  containing  five 
quantum  wells  of  varying  width  under  a  -2.0  V  reverse 
bias. 
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D.  Quantum  Well  Optoelectronics 
There  exist  a  number  of  applications  for  III-V 
heterostructures  in  the  area  of  optoelectronics  including: 
electro-optic  modulators  (2) ,  superlattice  avalanche 
photodiodes  (3) ,  quantum-well  envelope  state  transition 
devices  (QWEST)  (8) ,  and  multiple  quantum-well  optical  in- 


terference 

filters  (4) .  Two 

of  the 

most  promising  tech- 

nologies 

are 

presented 

here. 

these 

are  the: 

heterostructure 

quantum 

well 

laser. 

and  self 

electrooptical  effect  device. 

1.  The  Heterostructure  Quantum  Well  Laser 

The  optical  and  carrier  confinement  available  from  a 
heterostructure  quantum  well  provides  an  ideal  system  for 
a  semiconductor  laser.  The  physical  construction  of  a 
such  a  device  is  shown  in  Fig.  1.  The  device  is  composed 
of  two  heterojunctions  forming  an  undoped  GaAs  quantum 
well  bounded  at  either  interface  by  n  or  p-doped  AlGaAs. 
When  a  forward  bias  is  applied,  electrons  and  holes  are 
injected  into  the  well  from  the  n  and  p  AlGaAs 
repectively.  Population  inversion  is  achieved  in  the 
junction  due  to  the  excess  number  of  electrons  in  the 
quantized  conduction  levels  and  the  absence  of  electrons 
in  the  valence  levels.  This  process  is  shown  in  Fig.  2. 


In  the  junction  region  electrons  can  recombine  with  holes 
by  falling  to  lower  energy  levels  and  emitting  photons. 
The  recombination  radiation  may  interact  with  valence  band 


i 


•  Electrons 
o  Holes 


Figure  1.  A  double  heterostructure  GaAs/AlGaAs  laser  (from 
Ref.  1). 


Figure  2.  a)  The  energy  bands  of  a  forward  biased 


GaAs/AlGaAs  laser,  b)  the  index  of  refraction  profile  and 
c)  the  intensity  profile  (from  Ref.  1) . 


electrons  and  be  absorbed  or  interact  with  conduction  band 
electrons  and  create  more  photons  by  stimulated  emission. 
The  laser  cavity  is  provided  by  the  cleaved  (110)  faces 
forming  the  mirror  surfaces  of  a  Fabry-Perot  etalon.  The 
index  of  refraction  and  intensity  profiles  are  also  illus¬ 
trated  in  Fig.  2. 

A  major  advantage  is  associated  with  the  quantum  well 
energy  quantization:  recombination  can  proceed  between 
electrons  at  a  fixed  energy  level  with  holes  at  a  fixed 
energy  level.  This  feature  is  shown  in  Fig.  3  where  the 
energy  as  a  function  of  the  density  of  states  for  elec¬ 
trons  and  holes  in  a  quantum  well  structure  are  superposed 
on  that  of  a  bulk  sample.  Shown  are  four  transitions  in¬ 
volving  the  step-like  density  of  states  that  are  sources 
of  possible  laser  emission  lines.  In  a  bulk  material,  the 
carriers  are  distributed  in  a  continuum  of  energy  levels. 
Thus  narrow  linewidths  are  not  as  realizable  with  a  bulk 
semiconductor  since  recombination  cannot  occur  at  fixed 


energies  (5) . 

Of  critical  importance  in  the  design  of  QW  lasers  is 
the  well  width.  For  d  <  100  A,  the  well  width  approaches 


1 


the  electron  scattering  path  length,  1^.  Under  this  condi¬ 
tion,  electrons  are  not  scattered  sufficiently  to 


thermalize  in  a  large  density  from  the  conduction  band 


f, 
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Figure  3.  Density  of  states  diagram  for  an  AlGaAs/GaAs 
quantum  well  (from  Ref.  6). 


Figure  4.  Schematic  of  the  quantum  well  SEED  (from  Ref. 
11)  . 


n-type  AlGaAs  layer  directly  to  the  p-type  AlGaAs  layer. 
The  electrons  then  recombine  with  the  holes  in  the  AlGaAs 


at  energies  higher  than  the  GaAs  transition  energies  and 

with  lower  efficiency.  By  growing  multiple  quantum  well 

structures  with  thin  barrier  layers,  this  limitation  can 

be  overcome  (6,7).  Quantum  coupling  (tunneling)  is  then 

present  between  neighboring  wells  which  makes  the  overall 

active  region  much  larger  than  1  .  It  is  less  probable 

P 

that  an  electron  or  hole  will  recombine  in  the  AlGaAs  bar¬ 
rier  region  .  Therefore,  if  an  electron  isn't  scattered  at 
one  quantum  well,  it  only  has  to  cross  a  thin  barrier 
layer  before  being  scattered  at  the  next  well.  Using  this 
technique,  excess  carriers  are  effectively  collected  in 
the  quantum  states  of  the  GaAs  wells  and  thermalized 
2.  Self  Electro-optical  Effect  Device  fSEED^ 

The  quest  for  optical  computers  capable  of  parallel 
processing,  switching  and  exceptionally  fast  operation  is 
receiving  much  attention  (8,9).  This  concept  couples  to¬ 
gether  the  areas  of  holography,  fiber  optics,  Fourier  op¬ 
tics  and  semiconductor  lasers,  and  envisions  the  utiliza¬ 
tion  of  the  superior  characteristics  of  photon computing 

(either  in  lieu  of,  or  in  addition  to  computing  using 
electrons) . 

To  achieve  bistable  operation  a  system  must  have  two 
stable  output  states  for  the  same  input  over  some  range  of 
input  values.  Two  methods  of  obtaining  optical 
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bistability  are:  nonlinear  optical  interference  filters 
and  nonlinear  positive  feedback  etalons.  The  Self 
Electro-optical  Effect  Device  (SEED)  is  an  optically 
bistable  switching  device  which  employs  the  nonlinearity 
of  a  MQW  modulator  and  photodetector.  The  operation  of 
the  SEED  is  illustrated  in  Fig.  4.  (11).  It  is  composed 
of  a  p-i-n  diode  containing  a  MQW,  a  series  resistor  (R) 
and  a  voltage  supply  (V^) .  To  provide  switching  capabil¬ 
ity,  the  incident  light  wavelength  is  chosen  to  be  near 
the  exciton  resonance  for  zero  voltage  across  the  diode 
(call  this  wavelength  A^^)  .  With  low  optical  power,  almost 
all  of  the  supply  voltage  is  across  the  diode  because 
there  is  little  photocurrent.  This  voltage  provides  an 
electric  field  which  is  perpendicular  to  the  interfaces  of 
the  MQW.  This  field  shifts  the  exciton  absorption 
resonance  to  higher  wavelengths  (the  Stark  effect)  and  the 
optical  absorption  becomes  relatively  low.  Increasing  the 
intensity  increases  the  photocurrent,  decreasing  the  volt¬ 
age  across  the  diode.  The  reduced  diode  voltage  moves  the 
exciton  resonance  back  toward  increasing  the  absorp¬ 
tion.  This  action  further  increases  the  photocurrent 
which  leads  to  regenerative  feedback  -  the  cycle  of  in¬ 
creasing  absorption  producing  more  photocurrent  leading  to 
decreased  voltage  -  and  switching. 

This  process  can  be  further  explained  by  viewing  the 
theoretical  implications  of  optical  bistability  with  re- 


of  a  me- 


spect  to  absorption  (10) .  The  output  power,  P^, 
dium  with  absorption.  A,  and  system  losses  (i.e.  reflec¬ 
tion  and  scattering) ,  p,  can  be  expressed  by 

Rr  =  f  Cl-A)P=  TP  <i) 

where  T  is  the  transmission  of  the  medium  and  P  is  the  in¬ 
put  light  power.  Also,  we  can  identify  an  excitation  pa¬ 
rameter,  N,  by  N  =  AP  which  indicates  the  total  amount  of 
power  absorbed.  In  this  case,  the  absorption  will  be 
nonlinearly  related  to  the  excitation.  Equation  (1)  is 
readily  plotted  against  N  for  various  values  of  P  as  shown 
by  the  straight  lines  in  Fig.  5a.  Included  in  Fig.  5a  is  a 
measured  plot  of  T  vs.  N  for  a  GaAs/AlGaAs  MQW  structure 
showing  the  nonlinear  relationship  between  transmission 
and  excitation.  The  intersections  labelled  A  through  D 
represent  solutions  for  particular  input  power  values. 
Lines  A  and  D  intersect  once,  lines  B  and  C  twice,  indi¬ 
cating  one  and  two  solutions  respectively.  The  transmis¬ 
sion  properties  can  be  used  to  determine  the  transmitted 
power  as  a  function  of  input  power,  shown  in  Fig.  5b.  The 
system  is  clearly  bistable  since  the  points  B  and  C  have 
two  different  output  states  for  the  same  input  power. 

A  SEED  has  been  designed  and  tested  at  room  tem¬ 
perature  using  a  MQW  consisting  of  50  layers  of  95  A  GaAs 

wells  and  98  A  undoped  AlGaAs  (11) .  The  device  exhibited 

.  .  .  .  2 
extremely  low  switching  energies  per  unit  area  (-4  fJ/um 

.  2 

incident  optical  and  -14  fJ/um  electrical) .  The  switch- 


MATERIAL  EXCITATION,  N 


A  B  CD 

INPUT  POWER,  P 


Figure  5.  a)  Graphical  solution  of  bistability  for  a 
GaAs/AlGaAs  MQW.  b)  Transmitted  power  vs.  input  power  for 
an  optically  bistable  device  (from  Ref.  10) . 
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ing  power  was  inversely  proportional  to  the  resistance  and 
switching  time  was  proportional  to  R.  The  power-delay 
product  was  approximately  1.5  nJ  for  each  value  of  resis¬ 
tance. 
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E.  Direct  Interband  Transitions  in  Semiconductors 
The  behavior  of  optical  spectra,  measured  by 
ellipsometry ,  is  governed  by  absorption  processes.  In 
this  section,  interband  transitions  for  direct  gap  semi¬ 
conductors  are  considered  by  deriving  the  expressions  for 
interband  transition  rate,  joint  density  of  states,  and 
critical  points.  These  concepts  are  then  applied  to  the 
GaAs  optical  spectrum. 

The  relations  for  interband  transition  effects  are 
derived  from  first  order  perturbation  theory  (1,2,3).  The 
perturbation  approach  is  used  since  an  incident  light  beam 
produces  an  electric  field  in  the  semiconductor.  The  ma¬ 
jority  of  this  section  is  taken  from  Wooten  (1) . 

The  electric  vector  potential  (A)  of  the  applied 
light  wave  can  be  expressed  as 

jfs  -L  A5^^exp[^L(^.r-a)t)]  +  exp[”U{*f (i) 

where  a^  is  the  unit  polarization  vector  and  q  is  the  wave 
vector.  The  Hamiltonian  of  a  semiconductor  system  under 
radiation  is  modified  to  include  a  perturbation  term 

H'(f,t)  (2) 

Here  H^  and  H'(r,t)  are  the  unperturbed  and  perturbed  com¬ 
ponents  respectively.  H  can  then  be  expressed  as  the  sum 
of  four  terms.  These  terms  represent  the  kinetic  energy, 
the  potential  energy  of  the  ith  electron  due  to  the  ionic 
cores,  the  Coulombic  interaction  of  the  electrons,  and  the 
interaction  energy  of  the  ionic  cores,  respectively. 


where  the  index,  i,  is  summed  over  the  N  electrons  per 
unit  volume  (the  second  term  in  the  brackets  is  from  the 
relationship  between  energy  and  momentum:  p  =  E/c) .  The 
kinetic  energy  portion  is  obtained  by  taking  the  sum  of 
the  momentum  contribution  of  each  electron  in  the 
unperturbed  system  ("pj^)  and  the  modification  of  its  momen¬ 
tum  caused  by  the  applied  vector  potential  (A) . 

The  only  perturbation  component  in  equation  3  is  con¬ 
tained  in  the  kinetic  energy  term  since  the  other  terms 
involve  only  the  residual  electron-ion,  electron-electron, 
and  ion-ion  interactions  in  the  absence  of  an  externally 
applied  field.  An  expression  for  H'(r,t)  can  be  realized 
by  looking  at  the  kinetic  energy  term  of  H.  When  the 
squaring  operation  is  performed  on  this  term,  the  term 
is  associated  with  Hq.  The  term  containing  jAj^  is  dropped 
since  we're  only  concerned  with  the  linear  response  of  the 
system  to  weak  electric  potential  fields.  Thus  the  per¬ 
turbed  and  unperturbed  portions  of  the  Hamiltonian  are  de¬ 
termined  as 
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The  commutator,  [p,A]  is  given  by 


P*A^A‘?  =  "‘L‘^  V*  A  =  0 


since  the  vector  potential  polarization  is  transverse  to 


the  direction  of  travel. 


Therefore,  the  perturbed 


Hamiltonian  can  be  written 


*Tf\C 


Substituting  equation  (1)  into  equation  (6)  yields 

®  (So- V) 


where  the  first  and  second  terms  are  the  expressions  for 


the  incident  and  reflected  waves  respectively. 


We  can  write  equations  for  electron  energies  in  the 


conduction  and  valence  bands  using  the  energy  vs  wave¬ 


length  diagram  in  Fig.  1.  Here  the  zero  of  energy  is  taken 


at  the  bottom  of  the  conduction  band.  For  energies  in  the 


conduction  band. 


and  for  the  valence  band. 


£AM')-  -fa- 


W''''  -3  2.mK 

wher  m  and  m.  are  the  effective  masses  for  electrons  and 
e  h 


holes  respectively. 


In  considering  transitions  from  a  state  i  in  the  va¬ 


lence  band  to  a  state  j  in  the  conduction  band,  the  states 


are  described  by  Bloch  functions.  The  Bloch  functions  are 


compass  the  periodicity  of  the  lattice  (4) .  They  are  writ- 


=  IL"'''’- ,10, 

In  order  to  obtain  the  transition  rate,  the  perturba¬ 
tion  elements  of  the  optical  matrix  are  determined.  This 


is  done  by  looking  at  the  expectation  value  of  the  pertur¬ 
bation  elements  with  respect  to  the  Bloch  states.  Since 
we're  only  interested  in  absorption,  the  second  term  in 
equation  (7)  is  omitted  which  produces  the  following  for 
the  perturbation  elements: 

(11) 

The  periodicity  of  the  functions  u^  and  u^  allows  equation 
(11)  to  be  expressed  as  the  sum  of  integrals  over  the  unit 
cells.  This  provides  the  following: 

This  relationship  is  shown  in  Fig.  2  where  the  position 
vector,  r,  is  expressed  as  the  vector  sum  of  (distance 
to  the  nth  unit  cell)  and  r' (distance  from  the  unit  cell 
origin  to  a  point  in  the  unit  cell) .  The  perturbation  ma- 
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^  !>  ceil 

The  exponential  argument  in  equation  (13)  varies  rapidly 

and  periodically.  The  matrix  elements  will  be  very 

small  unless  (5) 

X'-X"+^  =  0  (14) 

which  also  serves  as  a  statement  of  momentum  conservation. 
Also,  the  wavelength  of  the  incident  light  is  much  larger 

0,  and 

(15) 

This  is  the  condition  for  direct  interband  transitions. 
Indirect  optical  transitions,  where  the  electron  changes 
momentum  during  the  optical  transition,  are  also  possible. 
The  change  in  momentum  involves  phonon  assistance  in  the 
form  of  either  phonon  absorption  or  phonon  emission.  The 
significance  of  equation  15  is  that  direct  optical  transi¬ 
tions  are  far  more  likely  to  occur  than  indirect.  By  sum¬ 
ming  over  all  the  unit  cells,  we  can  obtain  the  number  of 


(16) 


than  the  lattice  vector  length.  Therefore,  q 

A" 


unit  cells  in  the  crystal,  N^: 


2  R„]  _ 


Taking  V  as  the  volume  of  a  unit  cell  and  n  from  equation 


(10)  as  the  crystal  volume,  N  =  n/V. 

w 


Dl 


then  becomes 


(17) 


Here  the  second  term  in  the  integral  is  zero  since  Bloch 
functions  are  orthonormal  (i.e.  an  electron  cannot  be  in 


the  conduction  and  valence  bands  simultaneously) . 


Using  the  matrix  elements  of  the  momentum  operator. 


V  Oceij 

and  placing  this  expression  in  the  rewritten  form  of  equa¬ 
tion  (17),  the  perturbation  matrix  elements  for  direct 
interband  transitions  become 

Vli;  -  f  dr  =  ^  ffi  (15) 

'  ZmfV  'Jceil  ' 

The  transition  rate  for  an  electron  from  state  to  V’j  is 

wp  f  mv.\' s  (-) 

To  determine  the  transition  rate  from  the  total  number  of 
conduction  and  valence  band  states,  we  must  integrate 
equation  (20)  over  all  k-space 


We  can  apply  the  following  general  property  of  the  S 


function  to  equation  (21) : 


Ux-x^) 


This  application  produces: 


(23) 

Substituting  this  relationship  into  equation  (21)  gives 
another  form  for  the  transition  rate 


'<y  = 


where 


represents  all  the  k  values  satisfying 


.  These  k  values  represent  a  surface  in 


VV3VS 


5 


k-space.  We  can  alternately  integrate  over  this  surface 
with  respect  to  a  surface  element  dS: 


An  expression  for  the  density  of  states,  p(£,),  can 


be  derived  using  the  following  arguments.  Consider  two 
constant  energy  surfaces  in  k-space,  one  with  energy  £  and 
the  other  with  energy  £  +  de  as  illustrated  in  Fig.  3. 
If  a  cylindrical  volume  element  of  surface  area  dS  and 
height  kj^  (k^  is  the  perpendicular  distance  in  k-space 
between  the  two  surfaces)  is  used. 


since  the  density  of  states  in  k-space  is  (45r^)”^  per  unit 
volume  of  crystal  including  spin.  Also, 


d  £  -  £  •  dl  A  “  1  ^  ^  I 


Therefore,  the  density  of  states  may  be  written 

p(i)  =  JL.  C  (28) 

This  expression  can  be  substituted  into  equation  (25)  to 


obtain  the  joint  density  of  states 


Here  E ^  is  the  energy  difference  between  two  electron 
energy  bands. 

The  points  in  k-space  where  j—^O  are  known  as 
critical  points  or  van  Hove  singularities.  Optical  spec¬ 
tra  will  usually  have  structure  at  energies  corresponding 
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to  these  points.  It  is  evident  from  equations  (24)  and 
(29)  that  the  transition  rate  and  hence  the  joint  density 
of  states  are  highest:  1)  where  the  £  ^  and  the  £  bands 
have  the  same  slope,  2)  when  both  £  ^  and  £ ^  have  a  local 
maxima  at  the  same  point  in  the  Brillouin  zone,  3)  when  £^ 
and  £ j  have  a  simultaneous  local  minima,  and  4)  when  a  lo¬ 
cal  maxima  and  minima  occur  at  the  same  point  in  different 
bands.  The  structure  may  not  necessarily  be  sharp  due  to 
background  from  other  allowed  transitions  in  the  Brillouin 
zone  with  the  same  energy. 

The  energy  band  structure  for  GaAs  is  shown  in  Fig. 
4.  It  is  readily  apparent  that  there  are  a  number  of  pos¬ 
sible  transitions  at  optical  energies.  Also,  moving  from 
the  direct  gap  position  at  the  r  point  towards  the  L  or  X 
points  there  are  regions  where  the  valence  and  conduction 
bands  are  essentially  parallel.  Both  the  fundamental 
bandgap  point,  and  these  quasi-parallel  band  regions  con¬ 
tribute  to  the  optical  spectrum. 

Fig„  5  provides  a  more  detailed  view  of  the  GaAs  en¬ 
ergy  bands  by  labelling  a  number  of  critical  points  asso¬ 
ciated  with  prominent  transition  regions.  These  points  can 
be  detected  with  modulation  spectroscopy.  One  such  tech¬ 
nique,  electroreflectance,  produces  a  first-order  de¬ 
rivative  of  either  the  energy  or  the  portion  of  the  di¬ 
electric  function.  The  differentiation  deletes  the 


constant  background  spectrum  and  provides  sharp  features 


corresponding  to  the  critical  points.  An 

17  -3 

electroreflectance  spectrum  for  n-type  GaAs  (10  cm  )  is 
illustrated  in  Fig.  6.  Transitions  occurring  in  the  r 
point  region  are  indicated  by  the 

E'q  +  Aq  critical  points.  Symmetry  points  moving  towards 
the  L  region  from  the  zone  center  are  denoted  E^()  and 
Ej^O  Critical  transitions  near  the  X  point  are  la¬ 

belled  and  E2  +  5* 
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F.  MODFET  and  Multiple  Quantum  Well  Energy  Bands 
The  energy  band  bending  and  subband  energy  levels  for 
both  MODFETs  and  multiple  quantum  wells  can  be  derived 
from  quantum  mechanical  and  classical  considerations. 
The  solutions  to  the  physical  picture  presented  by  these 
two  cases  are  complex.  Therefore,  approximations  are  made 
in  each  case.  These  approximations  maintain  the  integrity 
of  the  physics  involved  and  yield  a  quantitative  solution 
to  the  problem.  It  should  be  noted  that  section  1.  deals 
only  with  the  energy  vs.  distance  relationships  (E  vs.  x) 
and  that  section  2.  involves  the  dispersion  relation  of 
energy  vs.  momentum  (E  vs.  k)  to  determine  the  allowed  en¬ 
ergy  levels  in  the  quantum  wells. 

1.  MODFET  Energy  Bands  fE  vs.  x) 

Several  types  of  theoretical  models  have  been  used  to 
incorporate  the  quantum  effects  of  the  two  dimensional 
electron  gas  into  device  modeling.  Yoshida  discusses  four 
types  of  numerical  models  which  have  been  developed  to  de¬ 
scribe  the  electron  distribution  with  respect  to  position 
(and  hence  the  energy  states)  at  the  MODFET 
heterointerface  (9) : 

1)  Classical  model  in  which  the  subband  energy  levels 
are  neglected  and  Boltzmann  statistics  are  used. 


2)  Classical  model  in  which  the  subband  energy  levels 
are  neglected  and  Fermi  statistics  are  used. 


8 


3)  Self-consistent  quantum  mechanical  model  including 


many  body  exchange  and  correlation  effects  with 


local  density-  function  approximation. 


4)  Quantum  mechanical  model  using  the  triangular  well 


approximation. 


The  most  accurate  method  involves  the  self-consistent  so¬ 


lution  of  Schrodinger's  equation  and  Poisson's  equation 


(method  3) .  All  four  methods  produce  results  which  differ 


by  less  than  10%  for  electron  sheet  concentrations  of  less 


12-2 

than  1.5  X  10  cm  (9).  In  this  section,  the  spatial 


relationship  for  MODFET  energy  bands  are  determined  using 


method  4.  This  method  provides  both  the  overall  band 


bending  and  subband  energies  at  the 


MODFET 


heterointerface  as  a  function  of  position.  The  approach 


shown  here  is  taken  mainly  from  Landau  and  Lifshitz  (4) . 


The  band  bending  at  an  AlGaAs/GaAs  interface  is  de¬ 


termined  through  Poisson's  equation  which  relates  the  po¬ 


tential  V(x)  to  the  space  charge  density  p(x)  in  the  fol¬ 


lowing  manner; 

I  T“  ~  (1.1) 

da 

where  is  the  permittivity  of  GaAs.  If  it  assumed  that: 


(1.1) 


1)  charge  neutrality  exists,  2)  the  doping  concentrations 


are  much  greater  than  the  intrinsic  carrier  concentrations 


3)  the  bulk  materials  have  zero  charge  density,  and  4)  the 


dopant  density  is  uniform,  then  p(x)  can  be  expressed  (1) 


P(t0  J 


where  the  first  term  is  due  to  the  bulk  GaAs  and  the  sec¬ 


ond  term  reflects  the  charge  density  in  the  i  quantized 
levels  is  the  wave  function  and  n  is  the  charge  density 

in  the  quantized  region) .  Typically,  the  quantity  is 

14  “3  18  “3 

negligible  (10  cm  )  in  comparison  with  (10  cm  ) . 

When  only  considering  the  bulk  GaAs  term,  equation  1.1  has 

as  it's  solution  (2)  .  .  .  v 

Vu)  =  -^  (1.3) 

where  is  the  extent  of  the  charge  carrier  density  in 

the  GaAs.  Equation  (1.3)  certifies  that  the  band  bending 
in  the  bulk  GaAs  is  parabolic. 

Poisson's  equation  must  be  solved  simultaneously  with 
the  Schrodinger  equation  in  order  to  determine  the  energy 
subband  structure  of  the  system 

+  =  E '/'CTd)  (1.4) 

Figure  1  illustrates  the  potential  energy  scheme  used 
for  the  MODFET.  Here  the  quasi-triangular  potential  well 
approximation  is  used.  This  approach  assumes  that  the 
electric  field  is  constant  in  the  region  of  the  interface. 
The  quantization  takes  place  in  a  small  region  near  the 
interface.  Thus  the  AlGaAs  potential  barrier  at  x=0  is 
considered  to  be  infinite  with  no  wave  function  penetra¬ 
tion  into  the  region  x  <  0.  V(x)  is  thus  written  (3,4) 

\  x>o  (1-5) 

where  F  is  the  electric  field  strength.  If  this  approach 
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is  not  used,  then  an  iterative,  numerical  method  needs  to 
be  used  to  calculate  the  subbands.  The  associated  boundary 
conditions  for  equation  1.5  are 


0  ^0  \ 
X) 


(1.6) 

In  order  to  simplify  the  Schrodinger  equation,  a  change  of 
variable  can  be  made 

■2!L  ^ 

(1.7) 


X  -  =  C£ 

where  c  is  an  arbitrary  constant.  Equation  4  can  be  re¬ 


written  in  the  following  manner 


Now 


(1.7)  can  be  expressed  as  only  a  function 


(1.8) 


of  f 


with  c  =  (1.^) 

^  can^also  be  written  (this  notation  will  ne  used  later) 


_  ^  J  ^  __  /  <1  \  1 _  ^  T  _ _  ^  _ _  ^  ^  ^  ^ 


Equation  (1.9)  belongs  to  a  group  of  differential  equa¬ 
tions  solvable  by  LaPlace's  method.  This  method  applies 


to  any  linear  equation  of  the  form 

=  0 

yf\  so 


(1.11) 

The  general  solution  for  (1.9)  is  expressed  as 

K  ^  7E.(J:)  cL-t  (1-12) 

where  K  is  a  normalization  constant.  C  is  a  complex  inte¬ 
gral  with  the  path  of  integration  taken  so  that  the  inte¬ 
gral  is  finite  and  non-zero.  Z(t)  is  given  by 

S  )  <iiz  (1.13) 


?<*>  -  i 

and  P(t)  and  Q(t)  are  ~ 


85 


?(t)=  2 


Tn»o 


TA 


and  Qt+)=  (1.14) 

WS6  ^ 

The  coefficients  of  (1.11)  are  readily  determined  from 
(1.9).  This  yields  a^  *  a^  =  =  b2  =  0,  a2  =  1  and  b^  = 

-1  which  provides  Q(t)  and  P(t) 

b,t°-  = 

=  PiCt)  =  a^-t'  =  t' 


Now  (1.12)  is  expressed  .» 

For  the  integral  (and  consequently  V’)  to  be  bounded,  the 
path  of  integration  must  be  chosen  so  that  the  integrand 
is  closed  at  both  ends.  By  placing  C  along  the  imaginary 
axis  (t  =  iu)  the  integrand  is  finite  for  all  f.  The  wave 
equation  thus  becomes  .  ^3 


L  can  be  written  as 

ipLS)--K  f «,)  *  {“Cos (.%M  +  -^>3  )  dlM 


(1.17) 


which  can  be  written  as 

... 

(1.18) 

since  the  integral  containing  the  sine  term  is  zero  (inte¬ 
gral  of  an  even  function  that  is  integrated  symmetrically 
about  the  origin).  *(f)  in  (1.18)  with  K  =  -2i/(jr)^'^^  is 
an  integral  representation  of  an  Airy  function  (the  Airy 
function  is  further  discussed  in  reference  10) . 

Substituting  the  second  boundary  condition  in  (1.6) 
into  (1.18)  yields  .  .  . 

from  which 


WVWiWV, 


^ is  now  given  a  function  energy.  The  zeros  of  the  Airy 
function  will  satisfy  the  condition  of  (1.20).  Therefore, 
E  has  a  discrete  spectrum.  The  roots  of  the  Airy  func¬ 
tion  are  denoted:  0  <  a.  <  a_  <  . . . .  <  a  (n  =  1,  2,  3 
. . . ) .  The  quantized  energies  can  thus  be  related  to  the 


roots  of  the  Airy  function  by 


=C 


(1.21) 


(1.22) 


The  wave  function  can  now  be  presented  as  a  function  of 
the  roots  and  X  ^ 

According  to  Ref  (3) ,  the  Airy  function  can  be  ap¬ 


proximated  by 


This  expression  is  zero  when  the  sine's  argument  is  (n  + 
l);r.  Setting  the  argument  in  (1.24)  equal  to  (n  +  1)t 


gives 


(1.25) 


The  discrete  energy  levels  art  .ow  known: 

Viewing  this  relationship,  it  appears  that  the  spacing  be¬ 
tween  energy  levels  would  be  most  affected  by  the  doping 
concentration  since  all  other  parameters  are  constant. 
Also,  the  higher  the  value  of  n,  the  closer  spaced  the 


subband  energy  levels  become  due  to  the  exponent  in  equa¬ 
tion  1.26  being  2/3  ordered. 


2.  Multiple  Quantum  Well  and  Superlattice  Energy  Bands 
fE  vs. 

There  are  at  least  three  main  approaches  for  the  de¬ 
termination  of  MQW  and  superlattice  energy  band  levels. 
These  models  produce  solutions  based  upon  the  solution  of 
Schrodinger^s  equation  for  the  MQW.  The  superlattice  can 
then  be  viewed  as  a  MQW  with  a  relatively  thin  potential 
energy  barrier.  The  three  approaches  are: 

1)  The  Envelope  Function  Approximation  (EFA) 

2)  The  Tight  Binding  Approximation  (TBA) 

3)  The  Kronig-Penney  Model 

A  complete  presentation  of  these  methods  would  be  quite 
voluminous.  General  comparisons  between  the  different  ap¬ 
proaches  are  made  in  this  section.  The  Kronig-Penney 
Model  is  derived  here  because  more  complex  models  are  not 
necessary  to  illustrate  the  properties  of  MQW  structures. 

In  the  EFA  the  wave  function  of  a  superlattice  or  a 
MQW  is  expressed  as  a  linear  combination  of  Brillouin 
zone-center  bulk  stf  tes  with  a  slowly  varying  envelope 
modification  along  the  superlattice  axis  (perpendicular  to 
the  heterointerface)  (5,7).  The  EFA  has  been  found  to  be 
in  overall  agreement  with  experiments  (6) .  In  these  ex¬ 
periments,  comparisons  between  theoretical  and  measured 


The  coefficients  of  the  wavefunction  are  determined  by 
diagonalizing  the  matrix  elements  of  the  Hamiltonian  ma¬ 
trix  until  the  eigenvalues  correspond  to  the  state 
originally  assumed. 

The  TEA  has  been  found  ro  be  in  good  agreement  with 
the  Kronig-Penney  model  for  the  heavy  hole  valence  band 
states  near  the  zone  center  (5) .  Here  the  valence  band 
heavy  hole  states  are  nearly  parabolic  (a  condition  asso¬ 
ciated  with  the  Kronig-Penney  model) .  Calculations  involv¬ 
ing  both  the  TEA  and  the  Kronig-Penney  model  disagree  by 
less  than  a  millielectron-volt  for  heavy  hole  energies 
from  the  valence  band  maximum  to  -0.26  eV  below  that 
point.  The  agreement  between  the  two  methods  fails  for 
light  hole  and  conduction  band  states  due  to  band 
non-parabolicity  and  the  fact  that  the  Kronig-Penney  model 
fails  to  account  for  band  mixing  (as  does  the  EFA) . 

The  Kronig-Penney  model  will  now  be  discussed.  The 
potential  energy  band  vs.  position  for  the  MQW  is  shown  in 
Fig.  2.  The  barrier  potential  is  and  the  barrier  and 
well  thicknesses  are  b  and  a  -  b  respectively.  Here  the 
band  bending  is  not  shown  since  it  considered  to  be  only  a 
small  part  of  the  total  band  structure.  In  the  physical 
sense,  the  AlGaAs  and  GaAs  layers  are  thin  enough  that  the 
band  bending  is  small  over  the  distance  concerned.  This 
is  unlike  the  MODFET  case  where  the  GaAs  is  unbounded  on 


In  this  case  enei.gies  with  E  <  are  considered  in 
order  to  determine  energy  quantization  within  the  well. 
The  approach  is  to  solve  the  Schrodinger  equation 
separately  in  the  V(x)  =  0  and  the  V(x)  =  regions  and 
apply  the  boundary  conditions  for  continuity.  The 
Schrodinger  equation  and  its  associated  solution  for  each 


region  are  given  by  (8) 

V  ^  0  t.'u)' 


%  U)  =  I'  ♦  B  •,  P  = 


(2.1) 

(2.2) 


- - -0  -  -^(Vo-6)W(Tt)=C> 

.  ot=  (2.4) 

'  .  •k^ 


The  first  two  boundary  conditions  are  established  with  re¬ 
spect  to  the  Schrodinger  equation 

1)  f.Co)  =4  and  2)^  =  ^ 

djt  X‘t> 

The  other  two  boundary  conditions  are  based  upon  the  peri- 


(2.5) 


•  -rr 


dx 


djc 


odicity  of  the  potential. 

^  1  '/’v  (-b)  =  e'*-*  Vo  'b;  « ) 

Applying  the  four  boundary  conditions 

1)  A+6«C  +  D 

2)  l)8(a-6)  =  ^ 

3)  Ce*'>4D«’*= 


4) 


C*e ip -  Be 


(2.6) 


(2.7) 

(2.8) 

(2.9) 

(2.10) 


These  four  simultaneous  equations  can  be  solved  by  setting 
the  determinants  of  the  coefficients  of  A,  B,  C  and  D  to 
zero.  This  produces  the  following  transcendental  equation 


C^Jk^^  »  Si(vK(!*f)Sin p(a-b)  t  CosKCotb)  Co^tpCa-b)!  (2 .  ii) 


since  k  is  a  real  constant,  the  right  hand  side  of  the 
above  equation  roust  lie  between  1  and  -1.  The  allowed  en¬ 
ergy  values  are  those  which  cause  the  right  hand  side  to 
be  within  this  range. 

The  above  relationship  can  be  roathematically  siropli- 
fied  and  roore  easily  represented  by  applying  further  ap- 
proxiroations  to  the  Kronig-Penney  roodel.  This  roodel  im¬ 
poses  the  restriction 

P=  (^)bV^  H  Cons+an-f  (2.12) 

If  a  delta  function  approximation  is  used,  the  barrier 
width,  b-»0  and  barrier  height,  V^— ^  ,  and  the  range  of 
quantized  energy  levels  are  easily  found.  Applying  these 
simplifying  conditions,  (V^  -  E) — sinhob— ♦otb  and 
coshab-^l.  In  this  case,  (2.11)  is  written 

(|5a)  +  Cos  C^q)  (2.13) 

A  plot  of  the  right  hand  side  of  (2.13)  as  a  function 
of  ffa  is  shown  in  Fig.  3.  It  is  readily  seen  that  with 
increasing  P  the  energy  gaps  between  the  allowed  energy 
ranges  become  smaller.  Thus  the  quantized  energy  levels 
become  more  discrete.  This  is  qualitatively  the  case  in 
the  superlattice  and  multiple  quantum  well.  For  a  par¬ 
ticular  AlGaAs  composition  the  value  of  in  (2.12)  will 
remain  constant,  therefore  P  can  be  increased  by  main¬ 
taining  the  well  width  constant  and  varying  the  barrier 
width.  Thus,  small  values  of  P  are  representative  of  the 
superlattice  case  (with  the  associated  "spreading”  of  the 
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Figure  3 .  Allowed  energy  ranges  for  a  multiple  quantum 
well  using  the  Kronig-Penney  model  with  a  varying  P  param¬ 
eter  (from  Ref.  5). 


discrete  energy  bands)  and  larger  values  of  P  are  in¬ 
dicative  of  the  multiple  quantum  well. 
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Chapter  3. 

VASE  Studies  of  6aAs/AlGaAs  Heterostructures 


A.  Preliminary  Conuments 

Esssential  to  this  are  an  explanations  concerning:  1) 
how  the  angle  of  incidence  is  selected,  and  2)  the  avail¬ 
ability  of  samples  for  analysis.  The  explanation  of  <t>  se¬ 
lection  is  presented  now  to  prevent  repetition  since  this 
procedure  is  done  prior  to  each  study.  Also,  it  is  impor¬ 
tant  to  consider  the  process  of  sample  procurement  and  how 
this  affects  the  research  effort. 

Selecting  the  Angle  of  Incidence 

The  sensitivity  of  ellipsometric  measurements  is  a 
strong  function  of  <f>,  with  maximum  sensitivity  occurring 
at  the  wavelength-dependent,  principal  (pseudo-Brewster) 
angle  (1,2).  At  the  principal  angle,  suppression  of  ei¬ 
ther  the  p  or  the  s  component  of  the  wave  reflected  from 
the  sample  occurs.  VASE  allows  measurement  at  several 
angles  close  to  the  principal  angle,  as  opposed  to  spec¬ 
troscopic  ellipsometry  (SE) ,  in  which  measurements  are 
made  at  a  fixed,  single  value  of 

The  principal  angle  can  be  estimated  prior  to  mea¬ 
surement  by  modeling  the  assumed  structure  for  the 
sample  and  generating  three-dimensional  plots  of  or  A 
vs.  wavelength  and  <f>.  The  3-D  plots  for  V*  and  A  for  the 
structure  used  in  the  interfacial  roughness  study  are 
shown  in  Figs.  1  and  2  respectively.  Two  particularly  in¬ 
accurate  measurement  regions  occur  at  ^  near  0*  and  A  near 
180*.  Reproducibility  of  measurements  in  these  regions  have 
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been  poor  with  A  varying  by  as  much  as  5  .  The  V*  plot  in 
this  case  reveals  that  there  is  a  region  of  the  spectrum 
from  6300  A  to  8500  A  where  0  is  close  to  0*.  Therefore, 
measurements  where  4>  is  less  than  7 3® would  not  be  advis¬ 
able.  In  the  A  plot  a  handedness  change  in  the  polariza¬ 
tion  ellipse  is  seen  for  <f>  >  76*  .  This  handedness  change 
causes  difficulties  in  the  data  fitting  procedure  and 
hence  values  of  <f>  >  76*  should  possibly  be  avoided. 

The  effect  on  the  0  and  A  spectrum  of  changes  in 
layer  thickness  and  composition  is  revealed  through  3-D 
sensitivity  plots.  These  plots  contain  the  difference  (in 
or  A)  between  a  sample  containing  reference  parameters 
and  a  sample  with  perturbed  parameters.  A  sensitivity  plot 
of  Si)  is  shown  in  Fig  3 .  In  this  case  the  roughness  con¬ 
tent  of  the  GaAs  layer  in  the  sample  shown  has  been 
changed  from  10  A  to  20  A.  The  sensitivity  of  i)  to  this 
change  is  determined  by  subtracting  the  20  h.  i)  spectrum 
from  the  10  A  0  spectrum.  From  the  figure  it  is  evident 
that  i)  would  be  sensitive  to  a  10  A  increase  in  the  GaAs 
layer  roughness  at  ^  =  75t  Thus  for  this  sample  4>  =  75*  ap¬ 
pears  to  be  an  appropriate  measurement  angle  when  consid¬ 
ering  roughness  sensitivity,  handedness  change,  and  accu¬ 
racy. 

Obtaining  Samples 

The  VASE  method  is  a  valuable  diagnostic  tool  for 
semiconductor  growth  processes.  Currently  the  RHEED  (re- 
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Figure  3.  A  3-D  plot  showing  the  sensitivity  of  \(i  to  a 
change  in  roughness  layer  thickness. 


flection  high  energy  electron  diffraction)  process  is  the 
most  popular  in-situ  growth  monitoring  technique  for  MBE 
grown  materials.  VASE  is  a  possible  replacement  for  RHEED 
because  VASE  can  provide  additional  information  such  as 
optical  constants  and  alloy  composition  as  well  as  detect 
monolayer  growth.  Increasing  interest  in  VASE  measurements 
is  manifesting  itself  in  inquiries  from  growth  facilities 
at  Perkin-Elmer  Corp.,  Honeywell  Corp. ,  University  of  Il¬ 
linois,  Varian  Corp.,  and  the  Avionics  Laboratory  at 
Wright-Patterson  AFB  to  name  a  few. 

Much  of  the  demand  for  VASE  resulted  from  profes¬ 
sional  publication  and  offers  of  collaborational  studies 
by  the  UNL  laboratory.  All  of  the  samples  used  in  this  ef¬ 
fort  were  obtained  from  Dr.  S.  A.  Alterovitz  at  NASA  Lewis 
Research  Center  who  procured  them  under  NASA  a  grant 
from  D.C.  Radulescu  and  L.F.  Eastman  at  Cornell  Univer¬ 
sity. 
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INTRODUCTION 


Crucial  to  the  development  of  heterostructure  ICs  and 
opto-electronic  devices  is  the  ability  to  grow  GaAs/AlGaAs 
layers  of  high  quality  and  uniformity.  Layer  thickness 
variations  across  a  wafer  will  have  a  profound  effect  on 
threshold  voltage  in  MODFET  IC  devices.  The  abruptness  of 
interfacial  regions  determines  the  suitability  of  the  wa¬ 
fer  for  the  production  of  semiconductor  lasers.  Therefore 
a  requirement  exists  for  the  determination  of  layer  homo¬ 
geneity  and  interfacial  properties.  In  this  study  VASE  is 
used  to  address  these  important  considerations. 

A  layered  sample  (originally  grown  at  Cornell  Uni¬ 
versity,  labeled  #1723)  produced  for  the  later  development 
of  MODFET  devices  was  studied.  The  epitaxial  layers  were 
grown  by  MBE  on  a  liquid  encapsulated  Czochralski  (LEC) 
semi-insulating  substrate.  The  high  quality  of  the  mate¬ 
rial  used  in  this  structure  is  demonstrated  by  the  77K 
.  5  2 

electron  mobility  of  1.7  X  10  cm  /  V  sec  with  a  2DEG 
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sheet  carrier  density  of  3  X  10  cm  .  The  sample  used 
for  this  study  had  2.5  by  1.0  cm  length  and  width  dimen¬ 
sions. 

The  structure  is  illustrated  in  Fig.  1.  A  1  urn 
undoped  GaAs  buffer  layer  was  grown  on  top  of  the  sub¬ 
strate  followed  by  a  200  A  undoped  Al^  3^^^  yAs  spacer 
layer  (part  of  thickness  t^^)  .  The  planar  junction  of  these 


two  layers  forms  a  MODFET  heterointerface.  Next,  a  300  A 
doped  layer  of  AIq  3^30  ^As  was  deposited  in  order  to 
provide  carrier  electrons  to  the  heteroj unction.  A  300  A 
n^  GaAs  cap  was  the  last  layer  grown. 

EXPERIMENTAL 

In  order  to  obtain  good  sensitivity  in  the 
ellipsometric  measurements,  p-wave  suppression  was  used. 
Thus  measurements  were  taken  at  the  principal  angle  where 
A  is  90*  over  most  of  the  spectrum.  A  He-Ne  laser  (6328  A) 
was  used  for  alignment  and  setting  the  angle  of  incidence 
to  obatin  A=  90*.  This  sample  was  measured  at  two  angles 
(75.5*  and  76.5*)  near  the  principal  angle.  The  thickness 
of  the  buffer  layer  was  not  measured  since  it  is  farther 
from  the  surface  than  the  penetration  depth  of  the  light. 
Also,  the  200  A  AlGaAs  spacer  layer  and  the  300  A  doped 
AlGaAs  layer  were  modeled  as  a  single,  500  A  layer  of 
AlGaAs  since  the  inversion  analysis  isn't  sensitive  to 
doping  parameters.  Analysis  of  other  samples  by  this  re¬ 
search  group  has  yielded  information  about  layer  doping 
through  the  Franz-Keldysh  effect  (1,2). 

The  sample  was  aligned  with  the  long  axis  perpen¬ 
dicular  to  the  plane  of  incidence.  The  beam  spot  on  the 
sample  is  ellipsoidal  with  approximately  a  1  mm  short 
axis  and  a  4  mm  long  axis.  Measurements  were  made  at 
three  positions  on  the  sample.  Spot  1  is  near  the  point 
that  was  originally  the  center  of  the  2-inch  wafer,  spot  2 


is  in  the  center  of  the  sample  and  spot  3  was  located  near 
the  rounded  edge  of  the  wafer.  No  chemical  treatment  or 
rinsing  of  the  sample  (in  order  to  remove  surface  oxides 
and  hydrocarbons)  was  performed  prior  to  measurement. 

The  spectroscopic  measurements  were  taken  in  the 
3500  to  8000  A  wavelength  range  in  50  A  steps.  The  near 
infrared  points  were  much  noisier  than  the  others,  as  both 
the  Xe  lamp  light  output  and  the  multialkali 
photomultiplier  efficiency  are  lower  in  this  range.  In  ad- 
dition  the  values  of  A  above  7400  A  are  rather  low  (  <  30 
)  .  This  means  that  <f)  is  not  close  to  the  principal  angle 
in  this  wavelength  range.  Therefore,  the  A  data  in  this 
range  are  not  very  sensitive  to  the  model  parameters  of 
interest. 

RESULTS 

The  experimental  results  of  the  ellipsometric  param¬ 
eters  V*  and  A  at  the  three  spots  are  given  in  Figs.  2a 
and  2b  respectively.  The  main  feature  of  Fig.  2  is  that 
the  results  for  spots  1  and  2  are  identical  over  the  whole 
wavelength  range,  and  that  the  difference  of  spot  3 
spectrum  is  small.  This  exhibits  the  high  degree  of  homo¬ 
geneity  across  this  sample. 

A  least  squares  analysis  was  tried,  using  both  V  and 
A  over  the  entire  3500  to  8000  A  range,  to  obtain  the 
thicknesses  and  x.  However,  the  program  did  not  converge 
to  a  meaningful  result  with  respect  to  MSE  or  model 
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parameters.  As  the  experimental  results  for  il>  and  A  above 
7400  A  (especially  A) ,  have  a  much  higher  error  than  the 
lower  wavelength  results,  the  and  A  fits  were  restricted 
to  the  range  3500  to  7300  A.  The  results  obtained  by  the 
inversion  in  this  case  are  presented  in  the  first  three 
rows  of  Table  I.  Plots  of  the  calculated  and  experimental 
0  and  A  versus  wavelength  are  shown  in  Fig.  3.  The  MSE  is 
larger  than  the  results  published  in  a  previous  research 
effort  using  GaAs/AlGaAs  heterostructures  (3) .  The  main 
reason  is  the  difference  in  structures.  Here  we  have  a 
partially  doped  AlGaAs  layer,  while  in  (3)  the  AlGaAs  was 
undoped.  As  a  result,  the  Franz-Keldysh  effect  is  more 
prominent  in  the  present  case  and  should  be  accounted  for 
explicitly  in  the  structure  model  in  order  to  obtain  the 
best  possible  MSE  and  data  fits.  For  simplicity,  these  ex¬ 
tra  parameters  were  not  included,  as  the  main  objective  is 
the  comparison  of  a  couple  of  spots,  not  an  absolute  char¬ 
acterization  of  the  material.  As  a  result  of  the  large 
MSE,  the  90  %  confidence  limits  are  rather  large.  For  ex¬ 
ample,  the  confidence  limits  for  spot  2  are:  +/-  83  A  for 
t^,  +/-  51  A  for  +/"  11  A  for  t^,  and  +/~  0.02  for  x. 

Another  factor  contributing  to  the  large  confidence  limits 
of  tj^  and  is  the  high  correlation  (  0.965)  between 

the  two  layers  which  was  determined  during  analysis.  It 
should  be  noted  that  the  bulk  values  of  the  dielectric 
constants  and  <2^^^  were  used  in  the  inversion 
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analysis. 

The  results  for  the  thickness  and  the  Al  concentra¬ 
tion  X  from  Table  I  show  remarkable  consistency  with  the 
RHEED  values  (nominal  values  in  Table  I)  obtained  during 
MBE  growth,  especially  for  the  AlGaAs  thickness.  Also  ob¬ 
servable  is  the  excellent  homogeneity  obtained,  with  none 
of  the  estimated  parameters  changing  by  even  1  %  of  their 
value  at  the  center  of  the  sample.  However,  the  large  MSE 
values  can  give  the  appearance  of  an  inaccurate  comparison 
between  the  three  spots.  Therefore,  another  inversion  was 
made  using  only  the  experimental  results.  As  the  0  pa¬ 
rameter  is  less  affected  by  the  Franz-Keldysh  effect  than 
A,  the  entire  experimental  range  3500  to  8000  A  was  used. 
The  results  of  this  inversion  are  given  in  rows  four  to 
six  of  Table  I.  A  comparison  of  the  experimental  and  cal¬ 
culated  ^  results  for  the  central  spot  is  illustrated  in 
Fig.  4a.  For  completeness,  A  versus  wavelength  was  was 
calculated  using  the  parameters  of  Table  I  (row  five)  and 
the  result  is  displayed  in  Fig.  4b.  There  is  a  negative  A 
range  in  Fig.  4b  (above  7600  A),  shown  as  A  near  360*.  The 
major  conclusion  from  Figs.  4a  and  4b  is  that  the  discrep¬ 
ancy  between  the  calculated  and  experimental  values  of 
both  0  and  A  is  rather  large  at  the  AlGaAs  bandgap,  i.e. 
where  the  Franz-Keldysh  effect  is  very  important.  This 
0-only  inversion  reinforces  the  correctness  of  the  thick¬ 
ness  and  composition  values  obtained  from  the  data 
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Table  I.  Results  of  the  Inversion  Analysis  for  Sample 
#1723  (all  wavelengths  and  thicknesses  are 
in  A) . 


Soot  # 

MSB 

Tvoe 

X  Ranae 

X 

^1 

% 

% 

MSB 

1 

ip,  A 

3500-7300 

0.280 

504 

328 

32.9 

1.089 

2 

ip,  A 

3500-7300 

0.279 

507 

328 

33.0 

0.988 

3 

ip,  A 

3500-7300 

0.277 

499 

325 

31.4 

1.282 

1 

3500-8000 

0.265 

506 

320 

36.9 

.0155 

2 

3500-8000 

0.264 

509 

321 

36.8 

.0156 

3 

3500-8000 

0.265 

502 

317 

34.4 

.0156 

Nominal 

iilUtii 

***** 

0.3 

500 

300 

25.0 

**** 

fits  where  i>  and  A  were  used  simultaneously.  The  correla¬ 
tion  between  the  various  parameters  is  now  negligible,  as 
it  does  not  exceed  0.9.  The  new  90%  confidence  limits  for 
the  central  spot  are:  +/~  59  A  for  t^^,  +/”  27  A  for  t^, 
+/“  8  A  for  t^,  and  +/-  0.046  for  x.  The  absolute  value 
of  the  parameters  is  now  different,  but  the  comparison 
among  the  three  spots  show  exactly  the  same  features  as 
the  V  and  A  inversions:  x  is  constant,  the  AlGaAs  thick¬ 
ness  in  the  center  is  highest,  and  the  Gaas  thickness  at 
spot  3  is  1  %  below  the  other  two  spots.  It  can  be  con¬ 
cluded  that  the  comparison  of  thicknesses  and  A1  concen¬ 
trations  among  the  spots  can  be  made  on  a  scale  that  is 
one  order  of  magnitude  better  than  the  90%  confidence  lim¬ 
its,  provided  the  structures  are  almost  identical.  This 
result  is  crucial  for  any  future  three-dimensional  plots 
of  thicknesses  and  A1  concentrations  at  various  spots  on 
the  same  wafer. 
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C.  Characterization  of  GaAs/AlGaAs  Superlattices 
INTRODUCTION 

The  main  purpose  of  this  study  was  to  determine  the 
ability  of  VASE  to  measure  the  physical  and  optical  prop¬ 
erties  of  GaAs/AlGaAs  superlattices.  These  superlattice 
samples  were  grown  by  Cornell  University  for  later  forma¬ 
tion  of  MODFET  integrated  circuits.  Along  with  the  deter¬ 
mination  of  layer  thicknesses  and  composition  the  resolu¬ 
tion  of  excitonic  transitions  was  a  primary  consideration. 

The  observation  of  excitonic  effects  in  superlattices 
at  room  temperature  greatly  enhances  the  suitability  of 
their  use  with  other  opto-electronic  devices  (1) .  VASE 
provides  important  information  about  these  devices,  spe¬ 
cifically  the  index  of  refraction  and  extinction  coeffi¬ 
cient  at  excitonic  transition  energies  in  superlattices 
(2) .  VASE  is  advantageous  in  that  both  the  real  and 
imaginary  parts  of  the  dielectric  function  and  the  refrac¬ 
tive  index  are  readily  measured  without  the  use  of  a 
Kramers-Kronig  analysis. 

For  the  two  samples  used  in  this  study,  measurements 
were  made  at  ^  =  75.5*  and  76.5*  (for  sample  #2207),  and 
at  ^  =  76*  (for  sample  #2352)  which  are  very  near  the 
principal  angles  over  most  of  the  spectral  range.  The 
data  for  these  experiments  were  taken  in  the  3500  to  8000 
A  spectral  range  in  increments  of  25  A. 


The  nominal  structure  of  the  two  measured  samples  is 
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Figure  1.  Superlattice  structure  for  samples  #2352, 
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Figure  2a.  VASE  model  for  sample  #2352. 
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shown  in  Fig.  1.  The  layers  were  grown  by  molecular  beam 
epitaxy  in  a  Varian  II  machine  on  an  undoped  liquid  en¬ 


capsulated  Czochralski  (LEC)  GaAs  substrate  (3) .  The 
substrate  was  rotated  at  7  rpm  during  deposition  with  the 
deposition  rates  and  the  A1  to  Ga  ratio  calibrated  using 
RHEED  intensity  oscillations  (4) .  The  superlattice  buffer 
provides  a  high  quality  barrier  at  the  quantum 
well/superlattice  interface.  The  superlattice  GaAs  wells 
are  nominally  20  A  thick  for  sample  #2207  and  30  A  thick 
for  sample  #2352.  The  75  A  AlGaAs  section  is  an  undoped 
spacer  layer  which  separates  mobile  carriers  within  the 
150  A  GaAs  well  from  ionized  impurities  in  the  325  A  doped 
AlGaAs  layer.  A  single  atomic  plane  of  Si  provides 
"delta”  doping  for  the  quantum  well,  and  the  400  A  GaAs 
layer  at  the  top  is  used  as  a  surface  cap. 

MODELING 

The  model  used  to  analyze  the  VASE  data  is  shown  in 
Fig.  2a.  The  optical  constants  of  bulk  materials  are 
normally  used  for  the  modeling  procedure.  However,  in 
superlattices  the  energy  subbands  along  with  the  added 
complexity  of  wave  function  overlap,  nullifies  the  use  of 
bulk  material  properties  (5) .  Because  the  optical  con¬ 
stants  of  real  superlattices  are  not  independently  known, 
the  superlattice  was  modeled  as  a  single  Alj^Gaj^_^As  layer 
of  unknown  thickness  and  composition.  A  native  oxide 
layer  was  also  incorporated  into  the  model  on  top  of  the 


GaAs  cap. 

RESULTS 

In  order  to  obtain  the  best  possible  fits  to  the 
experimental  data,  a  number  of  modeling  approaches  were 
utilized.  The  results  of  the  modeling  sequence  for  both 
samples  are  shown  in  Table  I.  The  initial  best  fit  was 
over  the  entire  (3500-8000  A)  spectral  range  where  the  MSE 
was  minimized  for  both  and  A.  The  MSE  was  extremely 
high  (>  5C0)  in  the  first  runs  because  of  handedness 
changes  in  the  polarization  ellipse  which  are  not  detected 
in  the  data  acquisition  program.  Once  corrections  for  the 
handedness  were  incorporated  into  the  experimental  data 
(between  7450  and  8000  A  for  #2352,  and  between  7000  and 
7450  A  for  #2207) ,  the  MSE  values  dropped  to  approximately 
25  for  both  samples.  Due  to  the  shorter  wavelength  light 
being  absorbed  closer  to  the  surface,  the  analysis  range 
was  decreased  to  3500-6800  A  in  order  to  focus  on  the  lay¬ 
ers  above  the  superlattice.  Results  from  the  3500-6800  A 
analysis  were  used  as  input  parameters  for  the  rp  only  data 
analysis  over  the  full  spectral  range.  This  ensured  that 
the  fitting  procedure  would  be  more  sensitive  to  the 
superlattice  parameters. 

The  final  model  for  sample  #2352  with  corresponding 
layer  thicknesses,  composition  and  MSE  is  shown  in  Fig. 
2b.  These  values  are  from  the  ^-only  type  of  data  fits 
shown  in  Table  I.  Although  the  confidence  limits  are 


smaller  when  minimizing  the  MSE  with  respect  to  both  V  and 
A,  the  MSE  was  reduced  by  approximately  three  orders  of 
magnitude  by  going  with  a  V-only  matching. 

Plots  of  the  experimental  and  calculated  VASE  data 
for  sample  #2352  are  shown  in  Figs.  3a  and  3b,  and  for 
sample  #2207  in  Figs  4a  and  4b.  These  plots  show  the 
best  fits  obtained  between  the  calculated  data  and  the 
experimental  data  for  the  different  regression  analyses. 
The  data  are  matched  particularly  well  considering  that 
there  are  seven  variables  with  varying  degrees  of  correla¬ 
tion  between  them.  The  good  data  fits  are  evidence  that 
replacing  the  superlattice  with  an  "effective”  Al^Ga^^^^^As 
layer  is  a  valid  approach  for  this  specific  case,  where 
the  AlGaAs  barriers  are  thick  when  compared  to  the  GaAs 
quantum  well  thicknesses.  The  handedness  changes  are 
evidenced  in  the  A  plots  by  transitions  passing  through 
180  and  are  not  real  structure. 

The  broad  peaks  in  the  V  data  at  5800,  6500,  7200 
and  7900  A  for  #2352,  and  6000,  6800  and  7500  A  for  #2207 
are  mainly  the  result  of  optical  interference  effects 
from  the  superlattice  region.  Spectral  features  at  the 
shorter  wavelengths  are  dominated  by  the  top  layer  of 
GaAs.  However,  the  spectral  details  at  all  wavelengths 
are  influenced  to  some  degree  by  each  of  the  layers;  in 
particular  the  surface  oxide.  This  is  shown  by  Figs.  5 
through  8  which  are  discussed  below. 


Table  I.  Results  of  Experimental  VASE  Data  Analysis  (all 
thicknesses  and  wavelengths  in  A) . 


Sample;  #2352 

Fit 


Type  A  Range 

!i 

V»,A  3500-8000 

33.2 

418 

405 

0.37 

141 

5,529 

0.36 

rl>,A  3500-6800 

26.3 

429 

415 

0.35 

157 

5,417 

0.34 

3500-8000 

25.9 

437 

405 

0.37 

141 

5,464 

0.32 

Nominal  Values 

25.0 

400 

400 

0.30 

150 

5,500 

0.30 

90%  Confidence 

Limits  (+/") 

13.0 

89 

204 

.068 

104 

259 

.072 

Sample; 

#2207 

Fit 

Type  A  Range 

h 

^3 

^4 

^5 

V>,A  3500-8000 

15.4 

329 

239 

0.34 

151 

5,134 

0.22 

Vi, A  3500-6800 

25.1 

276 

303 

0.32 

139 

5,089 

0.17 

3500-8000 

24.7 

286 

275 

0.35 

153 

5,115 

0.21 

Nominal  Values 

25.0 

300 

400 

0.30 

150 

5,520 

0.30 

90%  Confidence 

Limits  (+/“) 

20.9 

134 

318 

.135 

240 

489 

.  191 

MSE 

24.9 

5.71 

067 

■klfklt 


MSE 

21.1 

1.37 

076 

hlfklt 


^!v: 
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Figure  2b.  Final  model  for  sample  #2352. 
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Figure  3a.  Psi  vs  wavelength  for  sample  #2352 
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Figure  4b.  Delta  vs  wavelength  for  sample  #2207. 


#2352 


with  varying  superlattice  thickness. 


Figures  5  through  8  were  made  by  fixing  the  thickness 
and  composition  values  found  from  the  analysis  of  sample 
#2352  and  sequentially  varying  one  parameter  at  a  time 
for  a  single  Figs.  5  and  6  exhibit  how  the 

superlattice  affects  spectra  in  the  6000  to  8000  A  range. 
It  is  apparent  from  Fig.  5.  that  increasing  the 
superlattice  thickness  causes  a  corresponding  increase  in 
the  amplitude  and  a  shift  to  higher  wavelengths  of  the 
spectra.  However,  an  increase  in  the  value  of 
superlattice  composition  is  seen  to  cause  an  increase  in 
amplitude  and  decrease  in  the  wavelength  shift  (Fig.  6.). 
Variation  of  the  thickness  of  the  upper  layer  of  GaAs  has 
a  pronounced  effect  on  0  between  4400  and  5200  A  (Fig. 

7 . ) ,  and  also  causes  reduction  in  the  peak  heights  of 
maxima  and  minima  at  higher  wavelengths.  Changing  the  ox¬ 
ide  thickness  uniformly  shifts  the  amplitude  of  i>  at  all 
wavelengths  below  7200  A  (Fig.  8.). 

The  growth  quality  of  the  super lattice  was  determined 
by  studying  interfacial  smoothness  between  layers.  In  or¬ 
der  to  evaluate  roughness,  a  Bruggeman  effective  medium 
approximation  was  performed  on  sample  #2352  with  a  20  A 
mixture  of  GaAs  and  Al^  3^30  ^As  in  the  model  between  the 
150  A  GaAs  layer  and  the  superlattice.  The  solution 

showed  no  appreciable  change  in  values  for  layer  thickness 

» 

or  composition.  This  resulted  in  an  increase  in  the  MSE 
of  0.3%,  suggesting  that  physically  there  is  not  a  mixture 
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Figure  6.  Generated  psi  vs  wavelength  for  sample  #2352 
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Generated  psi  vs  wavelength  for  sample  #2352 
with  varying  thickness  of  upper  GaAs  layer. 
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Figure  9.  Excitonic  spectra  for  sample  #2207. 
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layer  present.  Another  modeling  trial  incorporated  an 
AljO^  layer  between  the  same  layers.  This  analysis  was 
done  to  determine  if  roughness  existed  due  to  the  burial 
of  oxygen  impurities  (which  are  A1  composition  dependent) 
in  the  GaAs  portion  of  the  GaAs/AlGaAs  interface  (6) .  The 
solution  gave  the  same  oxide  layer  thickness  to  within  2  A 
with  no  change  in  MSE.  Thus  the  superlattice  growth  qual¬ 
ity  is  considered  to  be  uniform  and  interfacially  abrupt 
since  incorporation  of  a  GaAs/AlGaAs  mixture  or  oxide  into 
the  regression  analysis  worsened  the  fits. 

The  exciton  peaks  are  shown  in  Fig.  9  for  sample 
#2207  measured  at  10  A  intervals.  The  narrow  excitonic 
structures  seen  at  6530  A, 7500  A  and  7700  A,  are  respec¬ 
tively  the  second  electron  to  heavy  hole,  e-hh(2) ,  first 
electron  to  light  hole,  e-lh(l),  and  first  electron  to 
heavy  hole,  e-hh(l) ,  transitions  in  the  superlattice  quan¬ 
tum  wells.  The  exciton  energies  are  in  qualitative 
agreement  with  values  predicted  by  subtracting  the  exciton 
binding  energy  (9.1  meV  in  single  30  A  GaAs  quantum  well) 
from  the  e-hh(l)  transition  energy  of  a  quantum  well  with 
a  30  A  well  width  and  a  wide  Al^  3^3^  ^As  barrier  (7,8). 
This  analysis  yields:  Eg_jjh(i)  “  1*582  eV  (7,840  A)  versus 
7,700  A  in  Fig.  9.  This  approximation  is  reasonable  since 
the  barrier  width  is  large  in  comparison  to  the  well 
width.  It  follows  that  the  e-lh  (1)  and  e-hh(2)  transi¬ 
tions  are  the  next  sharp  structures  that  appear  at  higher 
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layer  thickness  values.  The  XTEM  and  VASE  results  are  con¬ 
tained  in  Table  II.  The  XTEM  values  compare  closely  to  the 
VASE  values  obtained  in  the  only  data  fit  (Table  I) , 
particularly  in  terms  of  the  overall  superlattice  thick¬ 
ness  (+/-  40  A) .  XTEM  was  unable  to  directly  differentiate 
between  the  native  oxide  layer  (t^^)  and  the  GaAs  cap  layer 
{^2^  '  Therefore,  the  layer  thicknesses  for  the  oxide  and 
cap  layer  are  entered  in  Table  II  as  a  single  value  (460 
A) .  The  VASE  analysis  used  the  XTEM  layer  thicknesses  as 
model  input  values.  The  alloy  compositions  of  the  AlGaAs 
and  effective  superlattice  layers  were  variable  since  XTEM 
doesn^t  determine  these  parameters.  The  superlattice 
thickness  was  made  variable  in  order  to  determine  the  ac¬ 
curacy  of  the  solution  for  the  two  compositions.  As 
evidenced  in  Table  II,  the  90%  confidence  limits  on  all 
three  variables  are  improved  when  the  XTEM  values  are  held 
fixed.  Also,  the  superlattice  thickness  and  the  two  compo¬ 
sitions  are  essentially  constant  with  only  a  small  in¬ 
crease  in  MSE  (from  0.067  to  0.083).  Thus,  in  going  from  a 
seven  to  a  three  variable  solution,  the  accuracy  of  the 
VASE  analysis  is  further  verified. 

CONCLUSION 

Superlattices  in  MODFET/superlattice  structures  can 
be  effectively  incorporated  in  the  VASE  modeling  sequence 
by  representing  the  superlattice  as  a  bulk  AlGaAs  layer  of 
unknown  alloy  composition.  Analysis  of  five  layers  and 


two  alloy  fractions  is  possible  if  data  are  taken  near  the 
wavelength-dependent  pseudo-Brewster  angle.  Room  tem¬ 
perature  excitonic  effects  are  observed  and  growth  quality 
determination  is  realized  using  the  VASE  measurement 
process.  Thus  the  analysis  of  these  structures  through 
VASE  can  provide  accurate  and  nondestructive  assistance 
in  the  production  of  semiconductor  structures  for 
optoelectronic  and  microwave  communication  devices. 


Table  II.  Comparison  of  XTEM  and  VASE  Results  for  sample 
#2352  (all  thicknesses  are  in  A) . 


Analysis 

Type 

^2 

l3 

^4 

^5 

MSE 

XTEM 

460 

420 

*** 

160 

5420 

*** 

*  * 

VASE: 

Model 

26 

434 

420 

0.37 

160 

5420 

0.32 

★  * 

Result 

*** 

•kitic 

*** 

0.36 

*** 

5430 

0.32 

.083 

90%  lim. 

(V-) 

Ifk-k 

*** 

0.048 

Itlfk 

164 

0.008 

•k  k 
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The  purpose  of  this  VASE  study  was  to  determine  the 
amount  of  interfacial  roughness  at  an  inverted  interface 
where  the  GaAs  layer  is  grown  on  top  of  an  AlGaAs  layer. 
The  roughness  problems  inherent  in  GaAs/AlGaAs  layered 
structures  have  been  eluded  to  in  the  portion  of  this  the¬ 
sis  concerning  MBE  growth. 

The  basis  of  the  study  was  to  show  if  a  thicker 
AlGaAs  layer  increased  the  thickness  of  the  interfacial 
roughness  layer.  The  previous  atcempt  at  interfacial 
roughness  determination,  made  during  the  superlattice 
characterization  study,  didn't  indicate  a  roughness  layer. 
The  lack  of  roughness  resolution  in  that  study  was  poten¬ 
tially  due  to:  1)  a  low  measurement  sensitivity  associated 
with  the  inverted  interface  being  buried  deep  in  the 
structure,  and  2)  the  large  amount  of  highly  correlated 
variables  in  the  regression  analysis. 

To  overcome  these  previous  difficulties,  fairly 
simple  structures  were  needed  and  fortunately  obtainable 
from  NASA  Lewis.  The  sample  NASA  #1A  had  the  following  top 
to  bottom  nominal  layer  structure:  50  A  GaAs  cap  layer/ 
300  A  AIq  ^  GaAs  buffer  +  GaAs  substrate 
while  NASA  #2A  had  the  identical  structure  except  the 
AlGaAs  thickness  was  600  A.  These  two  samples  were  thus 
seemingly  ideal  for  our  purpose  since  the  inverted  inter- 


face  was  within  50  A  of  the  surface,  all  layers  were 
undoped  (no  Franz-Keldysh  shifts  with  doping) ,  and  the 
number  of  measurement  variables  were  low  (one  AlGaAs  com¬ 
position  and  four  layer  thicknesses  including  the  rough¬ 
ness  layer  added  during  modeling) . 

EXPERIMENTAL 

Before  the  actual  VASE  measurements  were  performed, 
it  was  important  to  determine  if  VASE  could  theoretically 
ascertain  roughness  on  the  order  of  a  few  monolayers  (5  A 
20  A) .  It  was  felt  that  if  VASE  wasn't  theoretically 
sensitive  enough  to  measure  interfacial  roughness  of  this 
size,  then  empirical  measurement  wouldn't  be  justified. 
Therefore,  3-D  generations  using  MAWEMA  were  performed. 
One  of  these  was  shown  previously  in  Fig.  3  of  this  chap¬ 
ter  (Preliminary  Comments  section).  This  plot  of  ^versus 
A  and  A  was  made  using  the  nominal  values  of  a  sample 
obtained  from  NASA  (NASA  #1A) .  The  sensitivity  of  at  least 
+/”  2®  in  r(>  for  a  10  A  roughness  perturbation  -  although 

small  -  was  considered  substantial  enough  to  measure  with 
the  VASE  system. 

Both  samples  were  measured  at  73.5  ,  74.5  ,  and  76.5 
at  20  A  intervals  from  3000  to  6000  A  and  10  A  intervals 
from  6010  to  8500  A.  The  10  and  20  A  wavelength  increments 
were  used  to  obtain  the  best  data  fits  during  analysis  and 
to  observe  any  fine  spectral  structure  that  might  be 


MODELING  and  RESULTS 

The  data  fitting  analysis  was  first  performed  with  no 
interfacial  roughness  layer  inserted  in  the  model  between 
the  GaAs  top  layer  and  the  AlGaAs  layer.  A  disturbing 
aspect  of  the  NASA  #1A  analysis  was  that  the  solution  for 
the  50  A  GaAs  cap  layer  gave  a  near  zero  value  (0.107  A) 
during  a  i()  only  fit  over  the  3000-8500  A  range  with  a  cor¬ 
responding  0.0288  MSE.  In  the  next  data  fit  for  #1A,  the 
GaAs  layer  thickness  was  removed  from  the  model.  The  same 
MSE,  AlGaAs  composition,  and  other  layer  thicknesses  were 
obtained.  The  experimental  V*  and  A  data  were  then  compared 
qualitatively  with  that  of  other  samples  containg  GaAs 
layers.  The  spectrum  of  sample  #1A  failed  to  contain 
structure  associated  with  the  and  Ej^+A^  fundamental 
GaAs  optical  transitions.  It  was  therefore  concluded  that 
NASA  #1A  had  no  GaAs  cap  layer.  Another  sample  from  the 
#1A  wafer  was  obtained  from  NASA  and  remeasured  to  pre¬ 
clude  the  possibility  of  the  thin  50  A  GaAs  layer  having 
been  etched  off  during  surface  cleaning  prior  to  the  ini¬ 
tial  VASE  measurements.  The  GaAs  layer  again  was  analyzed 
as  being  nonexistent.  At  this  point  it  was  concluded  that 
the  GaAs  layer  had  been  omitted  during  MBE  growth.  The 
lack  of  the  GaAs  layer  meant  that  an  inverted  interface 
didn't  exist  and  that  a  comparative  study  between  the  two 
samples  could  not  be  performed. 

Despite  the  lack  of  success  with  sample  #1A,  sample 
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#2A  was  still  analyzed  to  ascertain  if  a  roughness  ex¬ 
isted.  Modeling  of  #2A  was  done  with  and  without  a  rough¬ 
ness  layer  over  a  3000  to  8000  A  range  (the  8000  -  8500 
range  was  dropped  since  it  was  noisy  and  structurally  un¬ 
interesting)  .  A  GaAs  layer  was  found  in  #2A  which  gave  a 
31  A  thickness  in  a  ^  only  modeling  sequence  .  A  10  A 
variable  thickness  containing  a  50%/50%  mixture  of  GaAs 
and  AIq  inserted  as  a  Bruggeman  EMA  roughness 
layer.  A  zero  value  for  the  roughness  layer  was  produced 
by  the  regression  analysis. 

CONCLUSIONS 

From  this  experiment  it  can  be  concluded  that  the 
following  are  possible:  1)  VASE  isn't  sensitive  enough  to 
measure  sub-monolayer  aggregate  layer  thicknesses,  2)  sur¬ 
face  effects  (i.e  Fermi  level  pinning)  could  dominate  in¬ 
terfacial  effects  when  the  interface  is  quite  close  to  the 
surface,  or  3)  there  is  no  roughness  layer  present  in  the 
samples  investigated. 


